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CHAPTER I 


INTRODUCTION 

(A) HISTORICAL SURVEY 
1 , Fluid Dynamics 

Inspite of present day/ emphasis on atomic and 
nuclear physics and tendency of the theoreticians to avoid 
concrete models/ mechanics still remain the corner stone of 
physics. The mechanics of fluids is the study in which the 
fundamental principles of general mechanics are applied to 
liquids and gases. These principles are those of the 
conservation of matter/ the conservation of energy and 
Newton's law of motion. By the use of these principles we 
are enabled not only to explain and bring into rotation 
observed phenomenon but also predict at least approximately/ 
the behaviour of the fluids under a set of specified 
conditions. 

Towards the end of the Ninteenth century the Science 
of fluid mechanics began to develop into two directions 
which had practically no points in common. On the one side/ 
these were the science of theoritical hydrodynamics which 
was evolved from Euler's equations of motion for a 
frict ionless / non-viscous fluid and which achieved a high 
degree of completeness. Since/ the results of this so-called 
classical sciences of hydro-dynamics stood in glaring 
contradicat ion to experimental results/ it had little 
practical importance. For this reason/ practical engineers 
developed their own highly empirical science hydraulics. The 
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science of hydraulics was based on a large number of 
experimental data and differed greatly in its methods and in 
its objects from the science of theoretical hydrodynamics. 

At the beginning of the present century L. Prandtl 
(1904) distinguished himself by showing how to unify these 
two divergent branches of Fluid dynamics. He achieved high 
degree of correlation between theory and experiment and 
paved the way for the remarkably successful development of 
fluid machanics which has taken place over the past 
seventyfive years. 

Fluid mechanics; one of the oldest branches of 
physics provides the f oundat ion the understanding of many 
other aspects of the applied sciences and engineering. It is 
a subject of independent interest; in almost all field of 
engineering as well as in astrophysics; biology; 

biomedicine; meterology; physical chemistry; plasma physics 

e 

and geophysics. The development of aeronautical; chemical 
and mechanical engineering. During the past few decadeS; on 
the one hand there was exploration of the space and on the 
other hand the study of fluid mechancis was stimulated and 
as a result it has become one of the most basic subject in 
engineering sciences at present. 

2 . Magnet ohydrodynamics (MHD) 

When a conductor carrying an electric current moves 
in a magnetic field it experiences a force tending to move 
it at right angle to the electric field. Conversely; when a 
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conductor moves in a magnetic field a current is induced in 
the conductor in a direction mutually right angle to both 
the field and the direction of motion. These two statements 
first enunciated by Faraday/ constitute the laws of 
electromagnetism. The first is principle of electric motor/ 
the second that of dynamo. This is nothing in them to 
suggest that the conductor must be solid. In fact 
suggestions have^ been made that the motion of the sea may 
produce perturbut ions in the earth's magnetic field. Further 
tidal waves sweeping up the estuary of a river will cut the 
terrestial lines of force and generate a current which can 
be detected in a cable connecting two electrodes placed in 
the river on opposite banks. 

In the case when conductor is either a liquid or a 
gas/ electromagnetic force will be generated which may be of 
the same order of magnitude as the hydrodynami cal and 
initial forces." Thus the equation of motion will have to 
take these electromagnetic forces into account as well as 
the other forces. The science which treats these phenomena 
is called magnetohydrodynamics (MHD)/ other variants of 
nomenclature are hydromagnetics magneto-fluid dynamics/ 
magneto-gas dynamics etc. 

Most liquids and gases are poor conductor of 
electricity. As a consequence their motions can normally be 
treated by principles of fluid-dynamics. However/ it is 


possible to make some gases very highly conducting by 
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ionizing them. For ionization to take effect/ the gas must 
be very hot at temperatures upwards of 5000 K or so. Such 
ionized gases are called plasmas. The material within a star 
is a plasma of very high conductivity and it exists within 
a strong magnetic field. Consequently we expert MHD effects 
to be realised in stars. Further at the engineering level/ 
experiments have been made for electric generation by 
passing the ionized gas between the poles of a strong 
electromagentic / so that an electric current would be 
generated at right angle to the magnetic field/ and so the 
direction of flow of the plasma/ the current being 
collected by two spaced electrodes at right to the direction 
of the current flow. At the present time MHD generators are 
not a practical possibility owing to the difficulties of 
producing suitable/ efficient and stable plasma/ and 
sufficiently retractory materials to withstand the high 
temperature of the plasma. 

In summary MHD phenomenon result from the mutual 
effect of a magnetic field and a conducting fluid flowing 
across it. Thus an electromagnetic force is produced in a 
fluid flowing across a transverse magnetic field and the 
resulting current and magnetic field combine to produce a 
force that resists the fluid's motion. The current also 
generates its own magnetic field which distorts the original 
magnetic field. An opposing or pumping force on the fluid 
can be produced by applying an electric field 
perpendicularly t6 the magnetic field. Disturbance in either 
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the magnetic field or the fluid can propagate in both to 
produce MHD waves/ as well as up stream and downstream-wake 
phenomenon. The science of magneto fluid dynamics is the 
detailed study of these phenomenon, which occur in nature 
and are produced in engineering divices. 

3 . Historical Survey 

The subject of MHD had itsorigin in the speculations 
of certain scientists concerning the magnetism of celestial 
bodies . 

In 1989 Bigelow guessed that there were magnetic 
bodies on the Sun. Hale and Babcocks later confirmed this. 
In 1918 Zarmov made the attractive suggestion that magnetic 
field of the Sun and other heavenly bodies may be due to 
dynamo action, whereby the conducting material of the star 
acted as the armature and rotator of a self-exciting 
dynamo . 

In the inter-war period the astrophysicists, L. 
Cowling and Ferraro began to explore the formal theory of 

e 

MHD and its impl icat ions , with other scientists, and 

engineers like Williams and Hartman performed simple 
experiments on the flow of the conducting liquids in the 
laboratory- 

In 1942 Alfven [1] discovered the simplest example of 
the coupling between the mechanical and magnetic field 
force in a highly conducting motion in an external magnetic 
field. He showed that this interaction would produce a new 
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kind of wave. He called it MHD wave. These waves run along 
the lines of forces with a velocity which varies directly as 
the intensity of magnetic field and inversely as the square 
root of the density of the liquid. Hence/ transverse waves 
can be propagated in an electrically conducting 
incompressible fluid pervaded by a magnetic field and energy 
can be transmitted on the large scale exchange of fluid 
elements as much as the case of a non-conducting 
incompressible fluid. 

To obseEve Alfven's wave motion in a laboratory 
Lundquist [25] used mercury as liquid which he placed in a 
stainless steel cylinder and applied external vertical 
magnetic field of about 1000 gauss parallel to its axis. Due 
to dissipation true standing wave could not be excited. In 
1952 [21] tried liquid sodium as medium. Chapman and Ferraro 
[5] in early thirties developed the theory of 'magnetic 
storms ' . 


Ever since the appearance of work of Alfven and the 
interest in the subject has been growing rapidly and 
consequently/ a lot of papers on MHD have appeared in the 
sequence of time. 

In recent years/ the study of MHD has gained and 
increased attention of mathematicians and engineers in view 
of wide applications in Plasma Physics and Astrophysics. 
Hartman [13] was the first to solve onthe flow of conducting 
fluid cross a magnetic field. He studied one-dimensional 
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flow between two infinite stationary insulated paralle 
plates and found the exact solution in closed form. The 
same work was extended by Hartman and Lazarous [14] who 
performed experiments with mercury. They took rectangular 
and circular pipe with transverse magnetic field and it was 
followed by an increase in pressure gradient. Nigam and 
Singh [35] have discussed that heat transfer for Hartman 
flow. Stability of such flows have been discussed by Lock 
[23]/ Michael [27]/ pavelock and Tersov [37] and Stuart 
[46]/ Shercliff [41] studied the steady MHD flow in pipes 
and obtained exact solution for the axial flow across in a 
rectangular pipe for high Hartman number/ and velocity 
profiles and found to degenerate into a core of uniform 
flows surrounded by boundary layers on each wall. Shercliff 
[42] solved the problem of circular pipe with conducting 
walls. He also calculated velocity profiles/ pressure 
gradient and induced potential difference. The work of 
Hartman and Shercliff was extended by Chang and Lundgren 
[4]. They took MHD ducts found in a straight-pipe of 
arbitrary cross-section and conducting wall in the presence 
of transverse magnetic field. 

An investigation of MHD axial flow in an annular 
region with radial magnetic field was taken by Globe [6] 
who found exact * solut ion in closed form . in the next paper 
Globe [7] analysed the effects of axial magnetic field on 
the transition from laminar to turbulent flow for a pipe 
using mercury. Kapur and Jain [17] discussed the 
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investigation of Globe to a great extent. 


As regards the study of unsteady flow of a conducting 
viscous incompressible fluid inthe presence of magnetic 
field, Lin [22] indicated a class of exact solutions of MHD 
equations for laminar flow while a theoretical formulation 
and investigation of unsteady flow was made by lady 
Zhenskaya and Solomika [55], Lud ford [24] and Regirer [38], 
Regirer considered the problem of unsteady MHD flow between 
to parallel walls with an existing transverse magnetic field 
neglecting the incudction effects outside the flow region. 


Gupta [9] studied the flow of an electrically 
conducting fluid part in an infinite porous flat plate inthe 
presence of uniform transverse magnetic field. Kakutani 
[16] re-examined the problem studied by Gupta [9]. But he 
found that : an asymptotic solution for velocity and 
magnetic field is possible only if there is suction at the 
plate and suction is strong enough to check the outward 
transfer of vorticity due to Alfven wave. In these 

investigation the 'Hall' effect were neglected. The Hall 
effects later on were discussed by Sato [39], Yamanishi [55] 
and Sherman and Sutton [44] , when the flow is through 
straight channel. Katagiri [19] investigated the Hall 

effects on magn«etic hydrodynamic boundary layer flow on a 

semi-infinite flat plate. 

Gold [8] has obtained an exact solution of steady 

one-dimensional flow of an incompressible, viscous, 
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electrically/ conducting fluid through a circular cylinder 
in the presence of an applied (transverse) uniform magnetic 
field. He obtained the solution for all values of Hartman 
number. The corresponding problem discussed by Shercliff 
[43] was approximated by Gupta [10]. However/ his solution 
was valid for small values of Hartman number. He has 
obtained the velocity profile/ an induced magnetic field 
magnetic field by using interaction procedure upto second 
order term in the Hartman number. Gupta and Singh [12] have 
extended the problem of Gold to the unsteady flow for equal 
values of Reynold and magnetic Reynold's number. Later on 
Pathak [36] solved the same problem for general values of 
Reynold and magnetic Reynold's number. Laminar unsteady flow 
in an annulus^ under radial magnetic has been discussed by 
Gupta and Mittal [11] whereas Reynold's and magnetic 
Reynold's numbers are equal. For the first time/ the effects 
of wall's porosity on pressure distribution and velocity 
were considered by Berman [2] who has studied two dimensional 
steady flow in a rectangular channel with porous walls and 
obtained a solution valid for small values of suction or 
injection. Later/ Sellars [40] and Yuan [54] extended the 
work of Berman for large suction and injection values 
respectively.^ Morduchow [32] discussed the laminar flow 
through a channel and circular tube with injection by using 
the method of average. Berman [3] showed that for suction 
through walls of a porous circular tube. There were two 
solutions for hydrodynami cal 1 y developed velocity profiles 
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for a certain range of wall Reynold's number and/ in other 
range there was no solution. Terrill and Thomas [48] more 
recently showed that there exists dual solution for both 
blowing and suction/ except in the range were Berman [3] has 
found solution. Verma and Bansal [50] have obtained an extra 
solution of Na V i er-St okes equation for the flow of the fluid 
between two parallel plates/ one at rest and another in 
uniform motion with uniform porous walls of different 
permeability for small and large suction. Terrill [47] gave 
a first order perturbation solution for the steady flow in 
a porous annulus of different permeability for small cross- 
flow Reynold's number. Huang [15] amended the above 
Terrill's paper by using the method of quasi -1 inerisat ion . 
Nanda [34] obtained exact solution of the Navier-Stokes 
equation and energy equation/ for the case of a steady 
state of flow of the fluid through the annulus with porous 
walls while the inner walls moving with a constant velocity 
parallel to axis when other is at rest. Kapur and Malik [18] 
and Sinha and Choudhary [45] have studied the steady state 
laminar flow of the fluid between the coaxial porous 
cylinders rotating with constant angular velocity. Verma 
[51] has studied the pulsating viscous flow superposed on 
steady laminar flow in an annulus with porous walls. Verma 
and Gaur [52] have studied the unsteady flow in a porous 
annulus, when flow takes place under the influence of a time 
varying pressure gradient from rest. Muhuri [32] has 
considered MHD flow between two parallel porous plates when 
one was given an impulsing or uniformly accelerated motion 



( 11 ) 


with suction at both the walls. Mathur (26) has 
considered the unsteady flow of viscous incompressible 
fluid between two porous plates under a transverse 
magnetic field/ when the rate of injection of the fluid 
at the lower plate is equal to the rate of suction at 
the upper plate. 


Highly special though the property of 
ir rotat ional i ty may seem to be/ it is given great 
practical importance by the consequence of Kelvin's 
circulation theorem that material elements of a uniform 
field set into motion from rest remain without rotation 
unless they move into a region where viscous forces are 
significant. Recently Mittal (28/ 29/ 39) has made an 
attempt to find or create the condition under which the 
MHD flow becomes or tends to become irrotational in 
channel and ducts of different cross-sections under the 
action of an applied external magnetic field. 


(B) Hydrodynamical Flows in Ducts 
(a) Definition 

The mechanics of fluids is the study in which 
the fundamental principles of general mechanics are 
applied to liquids and gases. These principles are those 
of the conservation of matter, the conservation of 
energy and Newton's laws of motion. By the use of these 
principles we are enabled not only to explain and bring 
into rotation observed phenomenon but also to predict 
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at least approximately/ the behaviour of fluids under a 
set of specified conditions. 

For mathematical simplicity in describing flow 
in which the influence of viscosity is small a 
hypothetical fluid having zero viscosity may be 
postulated. This is known as an ideal fluid. 


To make headway in the study of fluid dynamics 
we must have a simplified model of the real fluid. As a 
model we take a continuous fluid and attribute to it at 
each point and instant a density/’^a pressure/ p/ a 
velocity v, a viscosity and a temperature T. 

If is constant we have the case of 
incompressible flow of a homogenous fluid and it is 
particularly * applicable to the motion of a liquid/ 
whose compressibility is negligible. 


The viscosity of real fluid never vanishes. 
Even if the viscosity is small/ as in case of air and 
water/ its effects are prominent near a boundary. For 
from rigid boundaries/ viscous forces can in many cases 
be neglected. In some problems the assumption/ the case 
of non-viscous flow/ gives a reasonably satisfactory 
description of the flow of a real fluid far from a 
boundary. In other problems/ however/ the discrepancy 
between the solution t^= 0 and that for very small may 


be quite striking. 



( 13 ) 


(t>) Steady or Stationary Flow 

In general the parameters such as velocity/ 
pressure and density/ which describe the behaviour of a 
fluid/ are not constant in a particular set of 
circumstances. They may vary from one point to another 
point/ or from* one instant of time to another/ or they 
may vary which both position and time. 

Steady flow is defined as that in which the 
various parameters at any point do not change with time. 
Flow in which changes with time do occur is termed 
unsteady or non-steady. Steady flow is simpler to 
analyse than unsteady flow. In practice absolutely 
steady flow is the exception rather than the rule/ but 
many problems may be studied effectively by assuming 
that the flow is steady/ for even though minor 
fluctuations of velocity and other quantities with time 
do in fact occur/ the average value of any quantity over 
a reasonable interval of time remains unchanged. 

A particular flow may appear steady to one 
observer but unsteady to another. This is because all 
movements are relative/ any motion of one body can be 
described only by reference to another body/ often by a 
set of coordinate axes. Since the examination of steady 
flow is usually much simpler than that of unsteady flow/ 
reference axes are chosen, where possbile, so that the 
flow with respect to them is steady. 



the 
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In general, fluid flow is three-dimensit iona 1 in 
sense that the flow parameters-velocity , pressure 
and so on vary in all the three coordinate directions. 
Considerable simplification in analysis may often be 
achieved, however, by neglecting the coordinate 
directions so that appreciable variation of the 
parameters occurs in only two directions, or even in 
only one. 


So-called one dimensional flow is that in which 
all the flow parameters may be expressed as functions 
of time and one space coordinate only. The single space 
coordinate is usually the distance measure along the 
central line (not necessarily straight) of some conduit 
in which the fluid is flowing. For instance, the flow in 
a pipe is frequently considered one-dimensional. In 
reality flow is never truly one-dimensional because 
viscosity causes the velocity to decrease to zero at 
boundaries . 

In two-dimensional flow the flow parameter are 
functions of time and two rectangular space coordinates 
(say X and y) only. There is no variation in the z- 
direction and therefore, the same stream line pattern 
could at any instant be found in all planes in the fluid 
perpendicular to the z-direction. 


Axially symmetric flow is three dimensional flow 


e 
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such that with a suitable choice of cylindrical polar 
coordinates R, 9, Z every physical variable is 
independent of the angle 0 or with a suitable choice of 
sphereical polar coordinates r, Q, 0 every physical 
variable is independent of the angle 0. Axisymmetric 
flow/ although not two-dimensional in the sense defined 
above/ may be analysed more simply with the use of two 
cylindrical coordinates (x and y). 


( c) Equation of Continuity 

The equation of continuity is really a 
mathematical statement of the principle of conservation 
of mass. This gives us an important formula which 
represents the conservation of mass and it is written 
in the form 


q + 


D 

Dt 


0 


( 1 . 1 ) 


( d) Euler's Equation of Motion 

Consider^ a small rectangular element of non- 
viscous fluids with sides of length x, y, z. The forces 
acting on it are the pressure on its several faces and 
possibility a body force F per unit mass due to 
gravitational or some other external forces. The 
acceleration produced is (^) where q is the velocity 
at some point in its interior. The Euler's equation of 
motion is 

F - i-,P = (1.2) 
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which in its alternative form becomes 
F--\7p = -V(q) + fox )x (1.3) 

( e) Bernoulli's Equation of Motion 

e 

The result of intergration of equation (1.2) is 
P 

— + — + qz = constant 

^ 2 

or, if we divide by g and put^g = W/ the specific weight 
of the fluid/ 

P 

— + — + Z = constant 

w 2g 

This equation is generally known as the equation of 
pressure in sieady flow. It was discovered by Danial 
Bernoulli and it is of fundamental importance for the 
entire theory of hydro-dynamics. 

( f ) The Navier-Stokes Equations 

It is clear that when viscous forces are acting 
in a liquid in motion/ Euler's equation of motion will 
no longer hold good and that it must be replaced by a 
modified equation in which the viscous forces are 
represented by additional terms and these additional 
terms are in facts : 

i •Vts,. 5) + 5 

where is the coefficient of kinematic viscosity. In 
place of Euler's equation we then have the Navier-Stokes 
equat i on 

E _ -A? P ^ q) +''V q (1.4) 
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(g) Reynold *s Number 

If we consider the motion of an incompressible/ 

constant viscosity fluid in a gravity field/ then the 

inertia forces . t , ^ r^ tji u j • 

ratio IS called Reynold's number and is 

viscous forces 

denoted by R. Reynold's number is important whenever 
viscous forces influence fluid motions. 

( h ) Some Solvable Problems in Viscous Flow 

There is no general solution to the Nevier-Stokes 
equations. Never-the-1 ess there are some special 
problems which can be solved. Few of them/ related to 
incompressible fluids/ are considered here. 

( i ) Steady Motion between Parallel Planes 

The region 0 4 z C h between the planes z = 0, 

z = h is filled up with incompressible viscous fluid 
(Fig. 1.1). The plane z = 0 is held at rest and the plane 
z = h moves with constant velocity V^. It is required to 
determine the nature of the flow when conditions are 
steady/ assuming that there is no slip between the 
fluid and either boundary/ neglecting body forces. 

Let p(x/ Yi z) be any point within the fluid. 
Then the velocity ^ at p will be of the form 

i 


q 




( 1 . 5 ) 
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The equation of continuity 
= 0 


vq .q = 0, gives 


B> V 

■BY 


and from (1.5) and (1.6) we infer that 

q = ^ ^ j 


( 1 . 6 ) 


(1.7) 


with nobody forces/ the Navier-Stocks vector equation 
may be taken in the form 


+ ( q.q )T7q +V^q 

c) q 

since the flow is steady = 0 . Also 

(q.\7)Vq = ^ v(z)j = 0 


( 1 . 8 ) 


and 


2 - , 

V q = V ( z) 


Hence (1.8) gives 


0 

- ( ^P i + ^P j 

k) + ^ V 
5z 


equating coefficients of the unit vectors/ 

If = 0 

(1.9) 

^ 

'iy 


(1.10) 

c)p ^ 

'B z 

6) 


(1.11) 

equations (l.\0) 

and (1.11) show that 

p(y) . Hence/ 

(1.10) 

becomes 




dp(y) 

dy 

— d^V(z) 

dz^ 


(1.12) 
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The L.H.S. of (1.12) is a function of y only, the R.H.S. 

is a function of z only. Hence each is a constant. As 

the fluid is moving in positive y-direction, the 

pressure p(y). should decrease as y increases. Hence — ^ 

dy 

0 and so we take 


dp(y) 

dy 


h- 


d^V ( z) 

^ 2 
dz 


( 1 . 12 ) 


The L.H.S. of (1.12) is a function is y only, the H.R.S. 
is a function of z only. Hence each is a constant. As 

the fluid is moving in positive y-direction, the 

dp 

pressure p(y) should decrease as y increases. Hence — 

dy 


0 and so we take 

dp(y) 

Hy = 


.2 , , 
d V ( z , 


az 


where P ^7 O, solving for v gives 

P 


V( Z) = A + Bz 


- <iT^> 


(1.13) 


when z = 0, V = 0 and when z = h, v = V^hence, we find 


v( z) = 


{■ 


< 2 


(1.14) 


(1.14) shows that the velocity profile between the 
plates is parabolic. 


{ ii) Steady Flow through Tube of Uniform Circular 

Cooss-Section (Poiseuille Flow) 

Fig. 1.2 illustrates the steady flow of an 
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inviscid incompressible fluid through a circular tube of 
radius. P is a pointin the fluid having cylindrical 
polar coordinates (R/ Qi z) , referred to the origin 0 on 
the axis of the tube which is taken as z-axis. We assume 
that there are nobody forces. Then continuity 
considerations applied to an annular shaped element of 
radius, R,R + 5" R of the fluid indicate that the fluid 
velocity is of the form 

q=w(R)k (1.15) 


Let us take the Navier-Stokes vector equation in 


form 


^ P + ( q p-x (vxp) (1.16) 

p 

> q 

For steady flow we have ~ 0 


The equation (1.16) becomes 


0 





2±. r, _L °—SL , 1 dp 


d(Rw' ) 

-m — 


equating coefficient s of unit vectors give 





^ P _ ^ d ( Rw ' ) 

3z ~ R dR 


(1.17) 

(1.18) 

(1.19) 


(1.17) , 


(1.18) 

dp( z) 
dz 


show that 


_ ^ d 

= dR 


p = p ( z) 
[RW (R) ] 


so that (1.10) becomes 

( 1 . 20 ) 


The L.H.S. of (1.20) is a function of z only; the R.H.S. 
is a function of R only. Hence, each is constant. As flow 
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is supposed to occur in positive direction/ we suppose 

(££) 4 0. Take each side of (1.20) to be - P, where P is 

dz 

a positive constant. Then 


d(Rw* ) 
dR 


-PR 


Or 


R 


dw 

clR, 


A - 


PR 

T7^ 


Hence 


dw 

dR 


R 


1 

2 


PR 

/u. 


and so 


w(R) 


= B + A log R - 


1 

4 


PR 


( 1 . 21 ) 


and w is finite on R = 0. Thus we require A=0, Also on 

2 

R = a/ w = 0 since there is no slip. Then = h (Pa /jx) . 


Hence 

w(R) = % (£ ) (A^ - R^) (1.22) 
The form (1.22) shows that the velocity profile is 
perholic/ i.e./ the plot of w against R from R=0 to a is 
of parabolic shape. 


The volume of the liquid 

anysection per unit that is : ^ 

a l\Pa 

Q = ;^w(R).2TrR. dR = 


B 


discharged over 


(1.23) 


If p denotes the pressure difference at two points on the 
axis cf the tube distance 1 apart/ then p = p/1. 
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(iii) Steady Flow between Concentric Rotating Cylinder 
Fig. 1.3 shows two concentric infinite cylinders 
of radii a,b (b y a) with viscous fluid in between. The 
inner cylinder is held at rest whilst the outer is 
rotated with constant angular velocity. Let (R, 9t z) be 

the cylindrical polar coordinates at point P in the fluid 

e 

referred to some point 0 on the axis of symmetry. Then 
when the motion is steady the velocity at P will be of 
the form. 


q = RW(R) 9 (1.24) 

where w(R) is the angular velocity of the liquid at 
radial distance R. The pressure at P will clearly be of 
the form 


P = P (R) 


(1.25) 


This time we t-ake the Nav ier-Stokes vector equation in 
the form 


- q X (Vx q) = - 

c>r ^ 2 ^ 


(1.26) 


for steady flow 



Equation (1*26) becomes 

- Rw^R = - i ^ R + 


(R 




(1.27) 


Equating the coefficient of unit vector/ gives 
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1 §P = 
? dR 



Rw 




From (1.29), we have, on multiplying by 


(R^ ^) 

dR dR 


= 0 


whence 

dw 

dR 



(1.28) 

(1.29) 


Taking w = 0, when R = a, w = , when R = b, we find that 
w(R) = b^ (1 - a^R"^) (b^ - a^)"^ (1.30) 

The fluid velpcity at radial distance R has magnitude 

2 2 2 -1 2 2 -1 
q = RW(R) = b (R - a R ) (b - a ) 

Thus the viscous stress at this location is 

b^(l + a^R~^)(b^ - a^)"^ 

Hence, the couple per unit length on the inner cylinder 
due to viscous drag has moment. 

(2[ra X ) X a X (2^b (b-a) ) = 4Trab f(b - a) 

This result forms the basis of measuring ^ for some 
liquid using a rotating visco-meter. 

( iv ) Steady Viscous Flow in Tubes of Uniform Cross- 

Section 


Before leaving the subject of solvable viscous 
flow problems, we treat the case of incompressible 
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unaccelerated flow through a tube of any uniform cross- 
section. We neglect body forces. Thus with (dg/dt) = 0/ 

F = 0/ the Navier-Stokes vector equation becomes 

^ p - \-T^p = 0 (1.31) 

Let us work with fixes coordinates axeS; Ox, Oy , Oz with 
Oz taken parallel to the flow so that q = w(x,y)k« 


Then (1.31) gives 

p 


( 1 . 32 ) 


Equating the coefficient of the unit vectors/ 


t)p 

c)X 

= 0 



( 1.33 ) 

Le 

= 0 



(1.34) 

a z 

= h 


,2 

3 w. 

+ N 2^ 

ciy 

(1.35) 


Equating (1.33) and (1.34) allow p = p( 2 ) so that (1.35) 
gives 

f ^ x / Y ) ] (1.36) 

^2 ^ X D y 

The L.H.S. of (1.36) is a function of z only^ the R.H.S. 
is a function of x / y only. Thus each of constant /say - P/ 
the minus sign being taken as we expect P to decrease as 
a increases- Then the problem reduces to the solution of 
the partial differential equation 
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5 w(x ,y) 


+ 


w(x ,y ) 


-p/^ 


(1.37) 


subject to w vanishing on the walls of the tube. 

( i ) A uniquencess Theorem 

To obtain solution of equation (1.37) it is 
desirable to establish a uniqueness theorem. A form which 

is little more general than required here is following 

,2 2 

If — 5“ + — ^ = f(x,y) at all points (x,y) of 
^x <5y 

a region S in the plane Ox, Oy bounded by a closed curve 
and if f is prescribed at each point (x,y) of S and at 
each point then any solution w = (x,y) satisfying these 

conditions is unique. 


Part C 

MAGNETOHYDRODYNAMIC FLOWS 
MHD is the science which deals with the motion 
of a highly conducting fluid in the presence of a 
magnetic field. The motion of the conducting fluid across 
the magnetic field and the action of the magnetic field 
on these current give rise to mechanical forces which 
modify the flow of the field. It is this coupling between 
the electro-magnetic and mechanical forces which 
characterises phenomenon. 

( a ) Fundamental Equations in MHD 

A general discussion of MHD involves the 
formulation and examination of the following equations. 
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(i) the equations of the electromagnetic field 
(Maxwell’s equations)/ 

(ii) the mechanical equations embodying the effect of 
the electromagnetic forces as well as other 
forces on the motion/ 

(iii) the equation of continutity/ 

(iv) the equation of heat transport. 

Also an equation of state and equations 
connecting the physical properties of the material with 
pressure and temperature must be considered. 

Maxwell's equations [20] involves the following 
vectors and their rates of change; the magnetic induction 
vector §/ the displacement vector 5 / the electric field 
intensity E ^nd the density of electric content J. 
Denoting di-electric constant by k/ magnetic permeability 
by jLL and using electromagnetic units throughout/ we 
have 

B = ^ 5 / 5 = ki ( 1.38) 

Maxwell's equations are the induction equation 

LI = - curl E (1.39) 

^t 

and Ampere's equation 

curl H = 4TrJ + (1.39) 

The following important subsidiary equations also hold 

= 0/ div 5 = 4Trq 


div B 


(1.41) 
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where q denotes the electric charge density. The first of 
these equations signifies that there are no free magnetic 
poles. The second equation is Poisson's equation. 

Differentiating the second of equation (1.38) 

with respect to time t, substituting for from 

equation (1.40) and using the vector identity div curl 

H = 0 we find 

div j = - ^ (1.42) 

0 L. 

This is the continuity equation for electric charge. 

Let denote the density of the fluid/ v the 
velocity of a fluid element/ p the hydrostatic pressure 
and F the body force per unit mass. When the 
electromagnetic terms are included/ the equations of 
motion of an electrically conducting liquid take the form 

+ (v. ) v] = - grad p +^J x H + F (1.43) 

The force F includes viscosity terms and any external 
forces such as gravity. 

In the case of a uniform incompressible fluid is 
constant. Then the equation of motion can be written as 

= - grad (— + v^) + tl (curl H) x H 

•bt ^ ^ 

(curl v) X V + F (1.44) 

Displacement currents have been neglected here/ so that 
curl H = 4115. This is useful from of the equation of 

motion . 
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The equation of continuity of the fluid, 
expressing the conservation of matter, can be expressed 
in either of the equivalent forms 

^ div ( V) = + div V 

The equation of state and the heat equation. In the case 
of perfect gas, the equation of state can be written as 

D rp 

P = — (1.46) 

m 

where T is the temperature, R is the universal gas 
constant and m is the molecular weight in terms of the 
hydrogen atom. For adiabatic changes the (p, )- relation 

i s 

(3 X. 

~ (P/ ■ 0, i.e. P/^ = constant (1.47) 

dt 

where r = -s the ratio of specific heats at 

constant pressu'"® and -onstant volume respectively. 

To these must be ndded the heat equation 

„ dS« dW P d 1 

^ ^ dt" ~ dt " “div grad T + - i (1.48) 

where S denotes the . specific entropy, W the specific 

internal energy, the thermal '■onductivity , and 0 the 

total dissipative function. 

( b ) The Magnetohydrodynamic Approximation 

Since the megnetohydrodynamic equations combine 
the full complexity of Maxwell's equations with the 
fluid dynamic equatio's, it is clear tiat they will be 
extremely difficult tf solve in their general form. 
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Accordingly some simplifications must be introduced if 
we are to proceed at all. Such approximations are 
generally used in engineering MHD as also in questions 
relating astrophysics. 

A number of approximations can be made which 
will be valid over a rather wide range of conditions. 

The first is to neglect the displacement current 
in Maxwell's equations, in the event that unsteady 
motions are of interest. 

Another useful simplification is to neglect the 
current flow due to the transport of excess charge as 
compared to the conduction current. 


The third is that we may neglect the 
electrostatic body force in the equation of motion. 

These three approximations can be shown to be 
valid, i.e. displacement currents can be neglected in 
Maxwell's equation: current flow due to charge transport 

i 

is quite small and electrostatic body force can be 
neglected inthe equation of motion. All three taken 
together can 'be considered to be magnetohydrodynamic 
approx i mat ions . 


Forms of Basic Equations for Reference 


For reference, the equations of motion are now 
written out in details. For incompressible fluid 
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(viscosity assumed quasi-constant)/ the equation in 
vector form is 


N V — — 2 — — _ 

[£;— +(v. )v];=p +1^. V- + E + JxB 

-d t ^ ' f e 

which is modified from of equation (1.44). Here the 

force F (per unit mass) of (1.43) is replaced by + 

(JxB). B / H, J and E have their usual meanings as 
given above. Actually ^E represents the electric force 
present when where exists net space charge/ JXB is 

the interaction force between the current and the 

magnetic field which becomes the predominant force in 
magnetofluid mechanics/ is the gravitational potential/ 
is viscosity as identified from / the magnetic 

permeabil it ity . 


Cartesian Coordinates 


( ;> u 


h U 

TjT 


ii' + uiJL., 
i y s 2 ' 


V ay ii: 



5 X 

E + J B 

ex y z 

- J B 
z y 


^ V 
^ X 

~ + V 

b V S V 5 

a^y' ^2 


= 

- 



^v 

^-7- + W -7- 

Dy ^ z 


- 


+ E + J 

e y 

B - J B 

Z X X z 

(£)w 

it 

u 

^ SW 
c) X 

S w 

+ V ^ + w 

s y 

S w ) 
i z 
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^ P 
1 >Z 

+ ^ 

^ J)x2 

^ 7 " 

- 


b 

d t 

+ E 

c z + 

'^X^y — '^y^X 
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Cylindrical Coordinates 


,dv. 


dv. 




dz 






2 «»«, 

- T ~5 SJ") - 

r r 


fr- Je'i^-Jz«e 


,^''a dv V ^■'e ^''e 

(tt* + V © + © IE” '*' ^ 7 aT” 

'dr r ~ 0 © zoz 


''r'^©> 


« . i i£ + u ( 
r d© 


a^v. 


^ 1 dv^ , 1,2 2 dv, 

H 4 . i -, - ^ ^ i:. ^ ^0 + ”"2 ~~ 

f' ■ I ^ r ^ ^ ^2 ^ r^ d© 

^ d©^ 


'^© . 

■7 * 




d__ 

d© 


£ 4- j 3 - J i3 

e^© 2 r r z 


-av^ ^ vg ^ 


dv .pv Vg £ 

+ Vj. ^ + V^ ^) 


r d© 

dv. 


d^v_ , dv , 

= - If r dF- -*■ ”■ 


r^ a© 


a^v. 


a 


* e^z -^r^© “ '^©^r 


Sperical Coordinates 
dv 
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If ^ -ff 12 


(sin Q ^~) 


1 2v ^ ^''q 

L ■'!>'■ ■ ' ■ '!> £ “ r _ 2 ^ COS 6 

e ^2 ^2 ae " ^ ''e 


4— 4] - ^ 

r^sln e 
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''e 


V 1 cos © 
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r sin 6 


r'^sin^ © 


2 

dvi 1 d Vi 

(=i" e SF^) " 7^8 


^'’r ^ 2 cos 6 £l 8 

r 2 sin 28 d+ 


r sin © 


r sin © 


^ ■*• e '^r^e “ '^©^r 
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( i ) Two dimensional flow 

Suppose of fluid moves in such a way that at 
any instant the flow pattern in a certain plane is the 
same as that in all other parallel planes within the 
fluid/ the flow is said to be two-dimensional. If we 
take any one of the parallel planes to the plane z=0/ 
then at any point in the fluid having Cartesian 
coordinates (x/y/z) all physical quantities (velocity/ 
pressures/ density/ etc. )associated with the fluid are 
independent of z. 

In the few cases of two-dimensional channel 
flow found in literature; the cross-sections considered 
are rectangular or circular. Although other cross- 
sections/ that iS/ elliptical/ parabolic/ etc./ might 
be solved in special cases, the rectangular and 
circular geometries have not been even austively 
studied. ^ 

( ii ) One dimensional MHD Poiseuille flow 

One of the simplest problems in 
magnetohydrodynami cs / called the Hartmann problem/ 
concerns the steady viscous laminar flow of an 
electrically conducting liquid between parallel planes 
in the presence of a magnetic field. The solution to 
this problem given us insight into the MHD generator/ 
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pump, flowmeter, and bearings, and it forms the basic 
method for treating all viscous flow devices. 

(iii) MHD Couette Flow 

If, in a parallel-wall channel, one wall is 
stationary and the other is moving, we have that is 
generally referred to as couette flow. The difference 
between Poisseulle's and Couette flow arises by virtue of 
different boundary conditions for in Couette flow the 
fluid velocity at the top plate is non-zero where as Lt is 

t 

zero in case of Poisseulle's flow. The top and bottom 
plates here are assumed to be insulators. 

( iv ) The rectangular channel 

The general problem of MHD flow in a rectangular 
channel having walls of finite, non-zero conducting is a 
very difficult one, involving the magnetic field solutions 
not only in the conducting fluid but in space outside the 
channel as well. The problem of the steady flow of an in 
compressible, viscous, electrically conducting fluid 
through a rectangular pipe in the presence of a transverse 
applied magnetic field was solved by Shercliff [43] while 
a problem of a channel with ideally conducting walls was 
considered by Ulfyand [49] and later by Chang and 
Lundren [ 4] . 

(v) 

In the study of MHD flow in tubes of rectangular 
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e 


cross-section the applied magnetic field was normal to two 
of the walls and parallel to the other walls. The 
analogous problem for the circular piple/ where the 
applied field is in the radial direction, has been treated 
by Pai [20], Although it is possible to obtain a radial 
magnetic field between two concentric iron pipes, it is 
more interesting to discuss the problem of MHD flow 
through a circular pipes in a uniform magnetic field. This 
latter problem is two dimensional, where as the former is 
one-dimensit ional . The latter problem is a practical 
interest because of the occurence of pipe flow through 
transverse magnetic fields in MHD pumps, generators and 
flow-meter. 

{ vi ) Flow between two coaxial rotating cylinders with 

a radial magnetic field . 

We now consider the steady flow of an 
incompressible conducting fluid between two concentric 
rotating cylinders composed of an insulating material. A 
Uniform radial magnetic field is applied throughout the 
flow region. The cylinders terminate at perfect electrodes 
which are connected through a load. The cylinders are long 
enough that axial variations are negligible, and since 
there is no preferred orientation for stable laminar flow, 
there can be no angular variations. 


(e) 


Development of Study of Flow Problem s 

The historical aspects of the development of MHD 
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problems have already been referred to in part (a) of the 
present chapter. Here it will be more beneficial to sketch 
a brief account of the contributions of certain important 
researchers in order to take the problem more clear and 
approachable. 

(I) Laminar flow of a conducting viscous liquid 

between parallel walls with a transverse magnetic 
field [13] 

The steady rectilinear flow of an electrically 

conducting liquid along a uniform channel under the action 

of a uniform transverse magnetic field gives rise to some 

of the few linear problems of magnetohydrodynamics. 

Hartmann has considered a one-dimensional problem in 

wiiich a liquid is constrained to flow horizontally 

between two infinite horizontal stationary plane walls. 

Taking rectangular axes OXYZ such that the boundary walls 

have equations z = ± L, and assume that the fluid 

velocity is given by = (v{z), 0,0) and the magnetic field 

by B = (o,o,B ) where B is constant (Fig.1'4), he found 
o o o 

V M Cosh M - Cosh (M /L) 

V = (M cosh M-Sinh nf (1.49) 

where M is Hartmann number defined by 

1/2 

M = (1.50) 


The form of the Hartmann velocity profiles [13] 
for various values of M are shown in Fig. (1.5). The main 

transverse field is to generate electric 


effect of the 
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currents which retard the fluid in the central regions and 
accelerate the fluid near the boundaries thus flattening 
the parabolic profiles experiences in the absence of the 
magnetic filed. When M 1 , that is when the viscous forces 


are large compared with the electromagnetic forces, or the 


channel is very narrow, equation (1.49) tends to the limit 


3v^(L^ - 2^) 



(1-51) 


when M >'7 1, that is for flows in wide channels or in the 
presence of strong magnetic fields, equation (1.51) tends 
to the limit * 

V = v^ [1-exp - M ( L- \z\ )/L] (1.52) 

Thus the velocity is approximately constant except in two 
boundary layers with thickness of order L/M. This result 
indicates that the velocity gradient near wall is of order 
M/ L and thus it is to be expected that a large magnetic 
field may given rise to instability of the laminar motion 
near the walls. This, however, is not found to be the case 
in experiments; in fact that the presence of a 
sufficiently strong magnetic field appears to stabilize 
such a motion and tends to inhibit the on set of turbulent 
mot ions. 



There are relatively few linear problems in 
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magnetohydrodynamics, the theory of the motion 
present of magnetic fields of fluids which 


in the 
conduct 


electrical currents according to Ohm's law. One such 
problem is that of steady/ Incompressible; rectilinear 
flow alowng a uniform channel under a constant transverse 
magnetic field. It is a case of considerable practical 
interest because of the utility of induction flow-meters; 
which rely on the generation of a measurable potential 
difference in the fluid in a direction perpendicular to 
the motion and to the magnetic field. However; laminar 
flow occurs more rapidly under a magnetic field since 
turbulence tends to be damped by eddy currents. 


Hartmann has solved the one-dimensional problem 
where the flow is between two parallel walls; the field 
being vartually infinite in directions perpendicular to 
the imposed transverse field. His solution is given by 
equation (1.49). The form of the velocity profiles; drawn 
by Sharcliff; for various values of M are shown in 
Fig. (1.6). The imposition of the transverse field 
generates current which related the fluid near the centre 
of the channel and acclerate the fluid near the walls; 
flattening the normally parabolic velocity profile until 
finally there is a core of uniform flow and an 
exponential boundary layer; whose thickness is of the 
order of L/M, at the walls. When M is large; boundary 
layers are found to occur on every wall. 
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(v) Laminar flow of an electrically conducting 
incompressible fluid in a circular pipe [41] 

Now we consider the simple case of a laminar 
steady flow of an electrically conducting incompressible 
fluid in a non-conducting circular pipe under an external 


radial magnetic 

field 

which 

is independent 

of 

the axial 

position 

of the 

pipe . 

We assume that 

there 

is 

no excess 

elect ri c 

charge 

< e 

= 0) , 

that the 

pipe 

is 

infinitely 


long, and then both the velocity and magnetic field are 
axially symmetrical, i.e. independent of the angular 

position and without tangential components. Furthemore, we 
assume that velocity has only an axial component which 

depends entirely on the radial distance. We use the 

cylindrical coordinates r, 0 and z. Where z = 0 arbitarily 
closes. In the present problem we have the following 
relations : 

= qe = 0' " ^z^^^ 

“© = ^z " ^z 

P = p(r,z) 

f ) = ^ ( tvv.) =0 (1.53) 

7)t 

where subscripts r,0,z refer to components along r,0, z 
directions respectively, f (r) means functions of r only, 
and f(r,z) means fuctions of r and z both. ( ) means any 

quantity of the unknown. 
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After the velocity distribution is known/ 
Hz and P can be obtained immediately by simple 
quadrature. One interesting point is that the variation 
of pressure in the z-direction is always linear no matter 
whether there is a magnetic field or not. Only with the 
external magnetic field. H = 0, the pressure varies with 
the radial distance/ while in ordinary laminar Poiseulle 
flow/ the pressure is constant for a given section. This 
point is also similar to that of turbulent flow. 


The 

been solved 
velocity dist 


similar problem for two dimensional case has 
by Lehnert . The general behaviour of the 
ribution of these two problems are the same. 


e 
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CHAPTER II 




CHAPTER-2 


I 

2.1 • SOME SELF - SUPERPOSABLE FLOWS 

In tl^is Chapter some self superposable flows in 
Cylindrical ducts have been discussed. This paper has 
been accepted for publication in ACTA Ciencia Indica. 

In literature, we find the flows mainly in ducts 
with circular and rectangular sections. Recently Mittal 
(5.6) and his Co-workers (7,8,9,10) have studied the 
flows in ducts with some non-costomary shape. In this 
chapter we have tried to find some flows in ducts of 
cylindrical shape. Phenomenon of superposability and 
self-superposability having been proved to be the most 
appalicable techniques in solving the hydrodynamic and 
magnetohydro-dynamic problems have been used here in 
bringing out very important conciusions about these 
flows. 

2 . 2 ON SOME SELF-SOPERPOSABLE MOTIONS IN CYLINDRICAL DUCTS 

ABSTRACT 

In this paper an attempt has been made to find 
some self-superposable fluid motions of incompressible 
fluid in cylindrical ducts. No stress has been given on 
boundary conditions and the solutions thus determined 
contain a set of constants. Pressure distribution of some 
of such ii flows has also been discussed. The curves 

nil 
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along which the vorticity of the flow become constant 
have also been attempted. Some irrotational flows and the 
surfaces along which the flows show the tendency of 
irrotationality has also been attempted. The aim of this 
paper is to introduce a method of solving the basic 

equations of fluid dynamics in cylindrical system of 

coordinaes by using the property of sel-superposability . 

INTRODUCTION 

It was shown by Ballash (2) that when a flow 
with velocity q satisfies the condition 

curl ( q X curl q ) = 0 (2.1). 

it becomes self-superposable . In the present paper the 
authors have attempted some fluid velocities for which 

( q X curl q )= p (say) can be respresented as the 
gradient of a scalar quantity. Attention has been 

focussed on incompressible fluids only. For such fluids 
some velocities have been determined which satisfy the 
condition in cylindrical system of coordinates and for 
which p can be represented as the gradient of a scalar 
quantity 9^ (say). The solutions thus found must be of 
self superposable type. Each will have a set of constants 
which can be determined by the boundary conditions. 

It has further been shown by Ballash (3) that if 
q^^ and q^ are self superposable flows and q”^ + are 

also self superposafble then q^^ and are mutually 
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superposable . By using this property some more self 
superposable flows have been determined. Attemps have 
also been made to find some curves along which the 
vorticity of the flow becomes constant and also the 
conditions of irrotationality . 

In this paper a method will be introduced to 
solve the basic equations of fluid dynamics in 

e 

cylindrical coordinates by using the property of self 
superposability . Mittal (5,6) Mittal et (7, 8, 9 & 

10) have recently solved the questions for different 
coordinates systems. 

FORMULATION OF PROBLEM 


Let q X curl q = p (2.2) 

For incompressible fluids, we have 

div q = 0 (2.3) 


Now if a solution of equation (2.3) be found in 
such a way that p can be represented as the gradient of a 
scalar quantity say , i.e.. 

p = grad 4 ^ (2.4) 

It will give a self-superposable flow. It has 
been shown by Agrawal (1) that such solutions will also 
satisfy the equations of motion for a steady flow. For 
determining a flow of liquid let us consider the flow in 
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cylindrical coordinates (r,&,z). If q^ t)® "the 
components of q at any point (r,&,z) in cylindrical 
coordinates (11), then in order to make equation (2.3) 
integrable we may consider the following cases 

CASE-I Let q^ = 0, In this case equation (2.3) will be 
satisfied by a solution. 


q = 0 

q=AU (r)aW (z) 

q = B Uv (r) V (9) 

r 

Where U (r), Uj (r) are integrable fuctions of r. V (©) 
and W (z) the integrable functions of 0 and Z 
respectively and A and B the constants. 


For this fluid velocity, it can be shown that p 
can be represented by the gradient of a scalar quantity 
given by 


0 


A“ W^/ . dr + A^ / UU^ dr 


,2 „2 r .,2 

.3 


> r 7 2 2 

- B" V"J ^ dr + B V 


( ^ 

r2 

+ B^ jw de - ABUUfW’ J vd© 


ABUUiV (wdz+A^u^ [ WW'dz (2.6) 


here U(r), (r), V (0), W(z) are represented by U, 

V and W respectively and U', U^, V and W' represent 

their differentials. 
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By choosing different suitable sets of values of 
U, , V and W me may get a number of self-superposable 
fluid velocities. One of such veloctiy can be obtained by 
taking . 


= U = r, V = W = Z, A = B 
the fluid velocity will become 

q = 0 
r 

qo = a r 2 

q = A r e = Ae 
z yx 

^ ^ -2n 3 2^ 2 2^ 2^„2 2„2, 

and 0 = A [ 1 r z +2 r +r +0 - r 0 1 


(2.7) 


( 2 . 8 ) 
(2.9) . 


If U, V, W are constants, then 


= 0 
= 

= D, 


( 2 . 10 ) 


and 0 


£l 

2r2 


( 2 . 11 ) 


CASE-II 


(i) 


When qQ =0, In this case, the self superposable 
flows may be 


q^. = Ai Vy ( 9-) W ( 2 ) 


^0 


0 

i 


( 2 . 12 ) 
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" ®1 ^2 ^2 


and 0 is given by 


? = - Aj. Vj Vj J ^ dr - v| | 


2 2 Uj dr 


- ®l'' 2 ( dr + b2 U| ( V 2 v\ dg 


- Aj Wj Wl J d^ + AjEi Wj^ j Vj^ Vj dg 


- AiB^W^Uj 


Vi V2 dg 


( Wj Wi dz 


- A 1 Bj, Vj Vj Ui 


W]_' dz 


+ J'^l 


(2.13) 


(ii) 


q = ©Z 

1 


=o 

% = A;l0 


( 2 . 14 .) 


0 = aJ 


2 72 


Af _ Aj Z 


+ a2 02 2^ 

2 r2 


. (2.15) 
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(iii) 


= c^/r 


qe = 0 


(2.16) 


and 0 = const. 


(2.17) 


Case-Ill 
may be 

(i) 


When q = 0, some self superposable flows 

2 


q = 


A2V3 (6) W2(Z) ^ 


q^ = ^2 U 3 (r) W 3 (Z) 


(2.18) 


and 0 


q. = 0 


= B? w 2 


>2 wg 




U 3 dr - bI 


- A2 B2 W2 W3 1— dr 


2 ”3 

3 . 

r 2 


J “3 "3 


dr 


a 2 wl 


V 




3 


‘^2®2^2^3^3 J 


+ Aj V§ [Wj wl 


V 3 dQ 
dz 


(ii) 


+ B 2 U 3 f W wl dz 


3 3 

^2 § 2 


(2.19) 


qo = 

q = 0 


( 2 . 20 ) 


and 0 = 1/3 A? r^Z^- a? r2z2 - A^ Z^ log r 
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+ 1/3 hi ^3 ^ ^ 


2 

r2 z2 (2.21) 


(iii) 


" *^ 2 /^ 


q = D 
^0 2 


(2.22) 


and 0 = 05 log r 


(2.23) 


Case - IV When q^. = = 0 a possible solution of 

equation (2.3) is given by 


% 

q„ 


0 

0 


^3 ^4 ^4 


J 


Case - V When q^ = q^ 

flow is given by 

= 


(2.24) 

0 , the self - superposable 


qg = 
<3, = 



(2.25) 


Case - VI When q_ = 0 = q„, the self - superposable 

' ■ 2 XT 


flow is given by 
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q = 0 

r 


^0 



(2.26) 


In all the above cases U , V , W 

n n n 

(n= 1,2, 3, 4 ) are integrable functions of r, 9- and z 

respectively and A,b C &D (n = l,2, .,..5) are 
constants which may be determined by boundary conditions. 


SUPERPOSABLE FLUID MOTION 

It has already been shown that the hydrodynamic 
flows given by equations (2.5) and (2.25) are self 
superposable . It can also easily be shown that if q^ and 
q 2 are the two flows given by equation (2.5) & (2.25) 
then P for q^ + q 2 can be represented by the gradient of 
a scalar quantity. Thus q^^ and q 2 will be mutually 
superposable and a flow 

q^ = A V (9) W (Z) .. 

r 

= B U (r) W (Z) y 

Q^. = C U (r) V (9) J (2.27) 

r 

is possible. The same flow can be determined by 
mutually superposing the flows (2.12) & (2.26), 
(2.18) and (2.24). 
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PRESSURE DISTRIBUTION 


It is interesting to note that fl is nothing but 
Bernoulli function given by (2.4). 

r2 


0 


= ( 


2g' 


+ h = p/g 


(2.28) 


Where q, g, h and P denote velocity, 
acceleration due to gravity, height above some horizental 
plane of reference and the pressure head. 

It is well known fact that for an 
incompressible fluid the pressure head P is given by 
(2.4) . 

p = P/p + Constant 

r- itq 

Where P is the pressure distribution. Also if 
the motion of the fluid be steady and slow then the value 
of h can be taken without much loss of generality as r 
for the flows (2.5), (2.8) and (2.10), 9 for the flows 
(2.12), (2.14), ( 2 . 16 )and Z for (2.18), (2.20) and (2.22). 

Thus for 

% 

the flow (2.20), the pressure distribution is 

p = K, Z + 9^_z^ + 9^ ...(2.29) 

r 

Similarly for the flow (2.22) taking C 2 = 
we have 

P = Kg log r + Kg (1/r^+l) + K^q 


(2.30) 
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Where ( n = 1 , 2 , 3, 4 , 5 , 6 , 7 , 8 , 9 & 10) are 
constants. Similarly the pressure distribution for the 
other flows can be determined. 


VORTICITY OF THE FLOW 


It was shown by Ballash [3] that for a 
self-superposable flow, vorticity is constant along its 
stream lines. It T is a unit. tangent along a stream line 
then 

T X q = 0 (2.31) 

by equation (2.16) and (2.31) it can be shown 

that 


T - ( 0,— ,y 1 ^ ' 1 1 1 1 I ■iiii- ■ — 

Jl/Cf + rVDf 

Hence the vorticity of the flow (2.16) is 
constant along the curve represented by equation (2.32). 

Similarly the curves of constant velocity can 
also be found for other flows. 


IRROTATIONM.ITY 


Vorticity for the flow (2.10) can be calculated 


as : - 


A 

Pi $2 + 

r 


. (2.32 


(2.33) 
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It is clear from equation (2.33) that the flow 
(2.10) is not rotational. 

For the flow (2.22). 

^ = 0 

Hence the flow (2.22) is irrotational throughout. 

Similar conclusions can be drawn for the other 
flows discussed earlier. 




e 
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CHAPTER - 3 


MAGNETOHYDROSTATIC CONFIGURATION : 

3.1 In this chapter a brief account of magnetohydro- 
static problems has been given and force-free magnetic 
fields are ^defined. Definitions and scope of the 
subject have been given and explained. This chapter 
contain one research paper which is my own contribution'. 
The paper is "On Some Magnetohydrostatic configuration 
in Cylindrical ducts. 


3.2 MAGNETOHYDROSTATIC PROBLEMS 

In general the action of the magnetic field on the 
electric currents flowing in a fluid will give rise to 
unbalanced mechanical forces which set up material 
motions. However/ in certain cases the magnetic forces 
either vanish or may be balanced by fluid pressures. So 
that the conducting fluid is in equilibrium. When a 
fluid is at rest, so that v = O, the basic equations are 
divH=0 ...(3.1) 




H 


1 ^ V H 




grad p -= pF + 


4^1 


(curl H) xH 


.. (3.2) 
. . (3.3) 


Here H denotes the magnetic field, h denotes the 
magnetic permeability, p denotes the pressure, p denotes 
the density and F denotes the body force per unit mass. 
The problem of magnetohydrostatics may be divided into 
two broad classes according as Lorentz forces does or 
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does not vanish/ of particular interest is the case in 
which the magnetic field exerts no force on the 
material . 


A VARIATION PROPERTY OF THE MAGNETIC FIELD 

The lines of force of a magnetic field in an 
infinitely conducting fluid move with the field. It is 
clear that no matter how the field lines are moved about 
by the fluid motion/ not all configurations can be 
distorted into each other. The problem of what 
properties make different hydromagnet ic systems 

topologically equivalent is very difficult. Here we give 
a simple variational deviation of force-free magnetic 

fields/ due to L, Woltzer [4]. 

Usually, the vector potential A of the magnetic 
field H is defined in such a way that div A = 0, A is then 
unique. For present purpose. Let us ignore the 

unnecessary restriction on A. Since 

H = curl A ...(3.4) 

and = curl (v x H) ••• (3.5) 


we may write 



X 


curl A 


... (3.6 


For this particular form of the vector potential 


at 


= curl A 


0 


... (3.7) 


For a region*^( enclosed by a stationary closed surface 
and writting/ 
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we obtain 

dl 

dt 


W j A. curl A dl, 


= -^1“. Curl A + A.curl^^JdT 


Integrating by parts and using (3.7) we obtain 
_ =2 (^-.Curl A)dT = 0 


. . (3.8) 


Integrating equation (3.8) we obtain 

.curl Adi 


I = Ji4 A., 


constant 


..(3.9) 


For all vector potentials which are constant ( 


on the boundary. It follows that 


i5L ! 


7) A 
^t 


= 0 ) 


A and curl A d( is 


an invariant of the motion. Consider now the stationary 
values of the magnetic energy M/ subject to the 
condition we have 

" -W a( 


I curl A)^d( 


(3.10) 


Let represent a Lagrangion Multiplier (the factor 

l/S*!;! convenience only). Then since I = constant/ 

M is stationary if and only if (M-'^^) I is stationary. 

8 h 


Thus 


s Jli. 


[(curl A) "/A^curl A]dT = 0 

for all possible small changes in A such that A =0 
on the boundary. Thus 

[2(curl 2>A).(curl A) 

— curlA-0<: A. curl ^ A] dT = 0 
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Integrating by parts, since ^ A = 0 on the boundary, we 

obtain 

2 SJSji A. (curl curl A -o^curlA d T = 0 
for arbitrary "Sa it follows that 
curl curl A =0(!^curl A 
curl H = 0(. H 


or 


d A = 0 


. . . (3.11) 


The result holds for all magnetic field such that 
c)'a 

= 0 on the boundary. A magnetic field 


or 


"c) t 


which satisfies equation (3.11) is known as a force-free 

•MP ^pag, 

r ^ ^ (curl H)xH vani- 

magnetic field, since the Lorentz force 5-=^^ 

shes. Force-free fields are particular solutions of the 

i 

magnetohydrostatic equation (3.3) in which the magnetic 
field does not disturb hydrostatic equilibrium. Equation 
(3.11) is the equation which governs beltrami fields in 
hydrodynamics. 

3.3 MAGNETOHYDRODYNAMIC STUDY STATES OF AN INVISCID 
INFINITELY CONDUCTING FLUID 

If V denotes the fluid velocity and p the 
density, the steady states consist of all solution of 
the equations. 

aw wa *<** 

p(curl v) X V - ^ (curl H) x H 

2 — 

= -grad p - p grade (^V ) +pF ...(3.12) 

curl (v X H) = 0 ...(3.13) 

Here, the body force F will contain no frictional 
forces. The general solution of these equations is 
difficult and we shall content courselves by pointing 
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out that if p is constant, the configuration in which 



4T1 f (3.14) 

is always a solution of these equations if F = 0 (no 

2 

external forces) and if p + ^pv = constant. Note that 
these equations also reproduce the relation between v 
and H in hydromagnet i c wave, 

Chandrasekhar ( 2 ) and others have developed in 
full the axially symmetric form of equations (3.12) and 
(3.13) using stream functions for both the velocity v 
and the magnetic field H. 

3,4 ON SOME-MAGNETOHYDROSTATIC CONFIGURATION IN 

CYLINDERICAL DUCTS 
ABSTRACT 

In the present paper it is proposed to study some 
magnetic fields with conservative Lorentz force in 
cylindrical coordinates, neglecting the effects of 
displacement current. It has been shown that self- 
superposable flow of an electrically conducting 
incompressible fluid permeated by a magnetic field with 
conservative Lorentz force may constitute a magnetohy- 
drostatic configuration under certain conditions. Some 
such possible configurations have been attempted for the 
conducting fluid flowering in cylindrical duct and 
pressure distributions of these configurations have also 


been determined. 
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INTRODUCTION 

The flows of conducting fluids have been studied 
mainly in ducts which circular and rectangular sections 
when the walls of the ducts are non-conducting and the 
transverse magnetic field usually applied. In the 
present paper we shall deal with certain problems 
connected with the flow of conducting fluids in ducts 
having more general cross-sections and shapes for 
difference coordinate system. Phenomenon of 
superposabil ity and self - superposabil ity have been 
proved to be most applicable technique in solving 
hydrodynamic and magenetohydrodynamic problems. Here we 
have used these techniques frequently in some important 
conclusions about MHD duct flow. Mittal et al [7/8] have 
discussed some of such problems in ducts of different 
shapes and cross sections. 

Lorerstz force plays many important roles in MHD 
flows. Various authors [1,4,7] have represented it in 
various coordinate systems. It is also proposed to 
discuss some magnetic fields in cylindrical coordinates 
for which Lorentz force is conservative. Some magnetic 
fields whose Lorentz force can be represented by the 
gradient of a scalar quantity have been determined. 

The basic assumption that displacement current 
field being small compared to the prevalent electric 
field can be discarded, have been taken into account in 
equations of the problem. The magnetohydrostatic 
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con f igurat i ons can be developed when the magnetic fields 
with conservative Lorentz force will act upon self- 
superposable flow of electrically conducting fluids. It 
is attempted to find out the conditions of such 
magnetohydrostatic configurations. Some such type of 
magnetohydrostatic configurations have been studied by 
Mittal [7], Mittal/ Thapaliyal and Agrawal [8], Mittal/ 
Thapaliyal and Salam [9] Mittal and Khan [10] for the 
ducts of different shapes. 

LORENTZ FORCE 


Let the Lorentz force be represented by 
L = - (H X Curl H) (3.15) 

Here H and ^ denote the magnetic field and 
permeability respectively. Now, suppose that 0 be the 
required scalar so that we can express L as 


L = - ^2!— grad Q^ (3.16) 

Mn 

Writing the various terms of equation (3.15) & 

(3.16) in cylindrical coordinates (10), we get 

(3.17) 

- V. ^3.18) 
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T YY?- ^ (-'^'^6)1 

H ^ De "-^Z. ' (3.19) 


Here L^) and (^j-' are components 

of L and H respectively at any point (r, e, z) in 
cylindrical coordinates. 


Also 


grad 0 


= If ; 


I V J- 1± + ?i 

' ^ •p. 5© 




(3.20) 


Where i^, i^, i^ are unit vetors along r, © & z directions 
respectively ; 


Equating the values of various terns from (3.15) & (3.16) 

we get 


^ f v\ T ^ _ L (.H )A ^ ^ t>^'^6)3 

■t ^ OT_ •^y,. H 

<^7- 


. (3.21) 


(3.22) 


(3.23) 


Values of when susbtituted 


CYic\K-V^.ol® 4 x) 


. (3.24) 




Also the magnetic field H should satisfy the equation 


Div. 


H=0 


on 


2>e- “i-T J 


D-2. 


(3.25) 


In order te make equation (3.24) integrable we may 
consider the following cases : 


CASE-I When H^=0, In this case equation (3.25) will be 


satisfied by the field. 

o 


Wk-- 




\ 






t. - 


K 


(3.26) 


When F(r), Fj^(r) are integrable functions of r, ^(G’) and 
Z(z) the integrable functions of 0 and z respectively and 
and p the constant. For this magnetic field, it can 
easily be shown that Lorentz Force L can be represented by 
the gradient of a scalar quantity given by 


f ^ J ^ fcH) (H) <Avi 

^ 5 M"(ej y lo-iaa - o(/b'?CM.)^,(.H) :z-'c7)Js'(cs;=^^ 

■Ar o(.^\V^ j2C2L^zVx)<=^^ (3.27) 


Now it is clear that magnetic field (3.26) must have 
Lorentz conservative force as it is represented by the 
gradient of a scalar quantity ^ given by equation (3.27) 

By choosing different suitable sets of values of F (r), 
Fj^(r), ^ (O')/ Z(z) we may get a number of magnetic fields 
with conservative Lorentz Force, one of such field can be 
obtained by taking 


F(r) 

Y(9) 

Z ( z ) 

P 


F^(r) 


= r 


o6 


J 


(3.28) 


the magnetic* field becomes 


H 


H, 


0 


H. 


^ © 


(3.29) 


= cxI .rz 
= 

K 

For this magnetic field represented by equation (3.29) it 
may readily be shown that the Lorentz Force can be 
represented by the gradient of a scalar quantity 


'■ 4 = <L 


^ " X 7?" -V JiY Jr Yr> 








■] 


(3.30) 
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Similarly one more magnetic field with conservative 
Lorentz Force is given by 


(3.31) 


For this magnetic field, is given by 


f - 


... (3.32) 


When 0, in this case the magnetic field with 


consevative Lorentz Force is given by 




. . . . (3 . 35 \ 


and 0 is given by 






^ r.. .U. 

^ K , <. jxw^-z. 

>3Z 




(3.34) 
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(ii) 


H 


Yl 


H 


e 


0 


H 


^ B 

>v ' 


and 4 = 


b^Z 


V 

7 - 


x^0 




1^ 


X'-b ■"V^ 

,^t> 1- 

— '^l_B-?' -V 


X^Y} 




'b K 


>1 7- 

°^\ B "Z 

2. n>- 


(3.35) 


CD'V' 


I ,7-1- X 'b 7_ T. 

^ ^ -V ^'3,^ __ ^ -z 


X- 


?~ 3— -i^ 

-V ^_Z 
XAl}- 


(3.36) 


( iii ) 


H 


= M 


H, 


H 


and 


f = 


0 

s> 

constant . 


(3.37) 

(3.38) 


CASE-I I I When H^-0, some magnetic field with conservative 


Lorentz Force are given by 


L^: 




"TL 




K 




e 




\d. 


o 


(3.39) 
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and 




(ii) 


^ o^j_ (.-Z- ) ’Z.^(-7' 5 ^'(P) 


o\<S 




A 


-V |b- 




(3.40) 


H 


H. 


H 




= 


Yv 

rz 


(3.41) 


and 

. 7^ *X i 7. V Y— *7* 0. A - I > 7^ /nTV 

Z n-. ~i — ri 

^o<X. & ^ T- "7- 

-T-^ -A ( 3 . 42 ) 

H := VZ 

^ le 

»0 = 

* =0 A ( 3.43 

and c|> = S^ ^ (3.44) 

In the above discussion all the nine magnetic 
fields in cylindrical coordinates given by equations 
(3.26), (3.29), (3.31), (3.33), (3.35), (3.37), (3.39), 

(3.41), (3.43) have conservative Lorentz Force and these 

are pressure balanced fields and it is a well known fact 
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that in such fields the linoc, p 

une lines of force and current are 

sittuated in the surface = Constt. 

SOME MAGNETO HYDROSTATIC CONFir.riRflTTnNQ 

If ^ denotes the fluid velocity of an 
electrically conducting imcompressible fluid permeated by 
a magnetic field H and f the density; the steady state 
consists of all the solutions of the equation given by 
(3.18) 

<* 

f cixha C. ^ 

-'f (3.45) 

^ -o (3.46) 


Here the body force F will contain no frictional 
forces. Let'V^^o^g, ^ t)e the components of ‘Vat any point 
(^' O'/ z) in cylindrical coordinates. These components has 
already been discussed by [3] for different cases as 

(i) 

<\l ~ o 

‘V - A U VO 

e ■ 


(3.47) 
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( ii ) 


^ ^ O 

A- 


- A K Z 
cv^ ^ A & 


( iii ) 


^J, 







(iv) 


( v) 


{ vi ) 


, A '-.' W , 

H 

K 

^ ^ , G Z 

K '- J: 


"\J 


n 


0 


A/ 

'=\j 

^2. 


z - ^\^e j 


D 




1- 



(vii) 


^ - 
^ ' 

^ - 
d ~ 






c 


(3.48) 


(3.49) 


(3.50) 


(3.51) 


(3.52) 


(3.53) 
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( viii ) 


( ix ) 


>v - i. - 


^ (5 ^ ( 

(3.54) 

^ 0 ) 




1 

( 

J 

(3.55) 


Where U, , V^, are the integrable functions 
of r; V, V^, V 2 / are functions of & and W, / W 2 » 
the functions of z. Also A, A^, A 2 / B, B 2 / C, C^, C 2 / 

D, , D 2 are constants which may be determined by 
boundary conditions. 


Now, let us consider the flow of an electrically 
conducting incompressible fluid with velocity given by 
equation (3.47) acted upon by the magnetic field given by 
equation (3.26). If 


F i’n}, -- 


. \ tK'l ^ 

- "A 


~L L^) 

- VJ L-L) 

> "tL 

0 } 

- \1 (_® ; ' 


cL ~ ' 

1 / — ^ 

f\ aw.cl 


/ 

1 




I 

\ qcr f \ 

(3.56) 

Then the flow 

(3.47) 

and 

magnetic field (3.26) 


will satisfy the equation (3.46). Also we have 
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(3.57) 


Wallen [13, 14] has shown that equation (3.57) 

is always a solution of equation (3.47), if 


V + ^ p ^ 


0 (no external solution) and if 
z- 


Constt . 


(3.58) 


Thus the flow of an electrically conducting 
incompressible fluid with velocity given by equation 
(3.47) acted upon by the magnetic field given by equation 
(3.26) constitute a magnetohydrostatic configuration if 
there are no external forces and then equation (3.56) is 


satisfied and pressure distribution is given by 

r X W -V- 1 

p = Constant 


r o’- w'- + 
L ~ ~~ 


Similarly we can show that the self -superposable 


flows and magnetic 

fields 

given 

by equation 

(3.48) 

and 

(3.29); (3.494 and 

(3.31) ; 

(3.50) 

adn (3.33); 

(3.51) 

and 

(3.35); (3.52) and 

(3.37) ; 

(3.53) 

and (3.39); 

(3.54) 

and 

(3.41) and (3.55) 

and 

(3.43) 

constitute 

some 

more 


magnetohydrostatic configurations in cylindrical ducts. 
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Magnetohydrostatic configurations thus found may 
be of either pressure balanced type or in very special 
cases force frtee type [ 5 ]. in all these, no secondary flow 
be possible because the magnetic forces would just balance 
the centrifugal effects [12]. 
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CHAPTER-4 


SOME MAGNETOHYDRODYNAMIC FLOWS 

4-1 In this chapter rotational and irrotational 

motions have been briefly studied. Some properties of 
vorticity and irrotational motion and their scope have 
also been given. This Chapter contains my own 

contribution in form of two research papers. These papers 
are : - 

1- On the Vorticity of MHD Unsteady Hele-Shaw flow 

of Non-Newtonian Fluid. 

2. On the Vorticity of an Elastico-Viscous Fluid 

near an accelerated plate. 

4.2 Plane Flow :- 


The movement 

of ideal fluids 

consists of 

two 

types differing 

from 

each other both 

physically 

and 

mathematically . 

They 

are called 

rotational 

and 


irrotational (or potential) flow. 

The types of flow will be derived by considering 
a fluid inplane flow and it a small area in the vicinity 
of the arbitrary point O (x,y) . The velocity at point O 
is V (X/y), point o', a small distance dr from O, has 
velocity/ v' = v + dv with components 

U ' = u + ^ dx +£u dy 

;ix Jiy 


(4.1) 
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V ' — V + ^ V dx + V dy 

<iy 

For convenience, let 


(4.2) 


a = ^_u 
Sx' 

b = ^_v 

c = 1 0 d + a u ) 

^ dx dy 


and “ i (^ V - ^ u ) 

2 Tx 1 y 

The equation (4.1) 
the following way 


(4.3) 

(4.4) 

and (4.2) can be written in 


U' = u+a dx + c dy - ^dy (4.5) 

v' = v+b dy + c dx - “^dx (4.6) 


To- visualize the meaning of the individual terms 
of this expression consider a rectangular fluid element 
(Fig. 4.1) whose one corner is the point O and whose 
diagonally opposite corner is O'. Clearly, u and v are 
the components of translational velocity, i.e., they 
indicate linear motion of all parts of the elements 
without change of shape (if a = b = c = "^ = O) . In view 
of (4.3) the terms adx and bdy express the change of 
velocity along dx and dy, respectively. However, since 
velocity is distance divided by time, a dx and b dy 
define the rate of change of length, the stretching rate 
of the edge of the element along the two principal 
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directions , 


= iz ' ^2 = ^ (4.7) 

Bx ciy 

so that, with (4.3) 

^ ^ = 2a (4.8) 

^ ^ ?)x by 


defines the change in the angle between the two edges at 
point O form its original • 90° . In other words, the terms 
containing a, b, and c in equation (4.5) and (4.6) define 
the deformation of the element. 

If the element is considered rigid for a moment, 

so that a = b = c = 0, then (4.3) requires that ^ = - ^_u 

bx, 'hV 

or 2 because of (4.7). This does not mean that 

is zero. In fact equations (4.4) and (4.7) show that in 

> 

this case. 

= h - ^2 = 2 V, 

That is, the element is subject to rotation around axis 
parallel to the z-axis through point O. It should be 
noted that the angleUj^ has the dimension 1/time, i.e., it 
is a change of agnle per unit time or an angular 
velocity, therefore. 

% =i3, = 1 (Sv - ^u) 

^ 1 ^ 

is the angular velocity with which the rigid element 
moves around the axis through point 0. 


e 
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Although this discussion has been confined to 
plane motion, it can be extended easily to three 
dimensions by analogous methods. As a result, the most 
general motion of fluid is made up of translation, 
deformation and rotation, in movement in space, the 
single relation (4.4) is replaced by three relations 
about mutually perpendicular axes and can beexpressed by 


t - 

1 

(c> W 

-civ ) 

(4.9) 


2 

^y 

Sz 

■S = 

1 

2 

(c) V 

-^u ) 
iy 

(4.10) 

yv = 

1 

(5 u 

- c>w ) 

(4.11) 


2 

5 z 

■Sic 


Fluid motions in which one or more of the terms 
b and is different from zero is called rotational 
motion. The vector 

+ i>t+ k ^ (4.12) 

is called the vorticity vefctor • ^ its 

components along the three coordinate axes x, y, z while 
i, j, k, are the unit vectors in the direction of these 
axes. In vector analysis the vector 2^that is, the vector 
whose components are 2*^, 2>i , 2tj is called the curl of v 
and written curl v, so that the rotational vector ^ and 
the velocity vector 'v are connected by the relation, 

= _1 curl V 
2 


(4.13) 
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ihe direction of vector is determined by the momentary 
location of the axis around which the element rotates/ 
and its magnitude is 



A flow system is called irrotational if the rotational 
vector disappears completely. According to (4.9), 
(4.10) and (4.11) this mean that 

b w - ^ V = 0 

'bY bz 

" b u - b w = 0 

'b z X 

h u - ^ = 0 

2>x ^y 

and the velocities u, v, w can be expressed as 
partial derivatives of a function 

j? (x, y , 2,t) 

or u = 

boo 

V = 

W f 

which can really be verified by substitution in equation 
(4.14). In vector notation, equation (4.15) becomes more 
simply 




V = grad 0 


(4.16) 
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Since the vectors sum of the three partial differentials 
is the gradient of 0 The function 0 is the velocity 
potential. In steady flow 0 is independent of time and a 
pure point function. 


By substituting the values for u, v, and w from 
(4.15) the equation of continuity for an incompressible 
fluid. 


o^u + 

one finds that, 


civ + "SiW =0 

in irrotational flow 






t# 



0 


(4.17) 


Here "v is the Laplace operator. For this reason the 
continuity equation (4.17) is often called the Laplace 
equation. It has been shown that in irrotational flow, 
the velocity v can be derived from a potential. There 
force such flow is frequently called potential flow. 


EQUATION OF MOTION FOR THE VORTICITY 

We shall now obtain an equation of motion for an 
incompressible, inviscid fluid as 

^u - u = grad w (4.18) 

> t 

where 

w ~Eti. |u]^ + v (4.19) 

P 2 


Taking the curl of 


equation (4.19) we obtain 


This 


is 


- curl 
the required 


(u ) 
equation 


= 0 

of motion 


for . 


MOTION OF VORTEX LINES 


(4.20) 


It is worthwhile to note that the induction 
euqation for an infinitely conducting fluid. 

^ = curl (U X H) 

and the hydrodynamic vorticity equation 

= Curl (u xS ) 

B t 

obtained when the external body forces per unit mass are 
conservative and similar. It follows that the vortex line 
move with the fluid (Helmholtz' theorem) and if the 
vorticity is initially zero, it is zero for all time. 

4 . 3 THE SOURCE OF VORTICITY IN MOTIONS GENERATED 

FROM REST . 


It can be established that changes in the flux 
of vorticity^ across a material surface element take place 
solely as a consequence of local diffusion of vorticity, 
by viscosity. 

Vorticity flux or circulation can not be created 
in the interior of the fluid but once there it is spread 
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by the action of viscosity. 

This raises the important question of the 
ultimate source of vorticity in motions which have been 
generated from rest in a fluid of liniform density. 
Initially the vorticity is every where zero, and the 
motion must remain wholly irrotational unless vorticity 
diffuses across the surface bounding the fluid. Real 
fluid motions which can be seen to possess vorticity over 
atleast part of the field are common so that we are lead 
to expect that some mechanism exists for the generation 
of vorticity at the boundary of the fluid. 

When the fluid is bounded wholly or partly by a 
solid, any remaining part of boundary being at infinity 
where the fluid is at rest, such a mechanism is provided 
by the no-slip condition. Mechanism for the generation of 
vorticity do exist at other types of boundary, such as a 
free surface at which the pressure is constant and the 
tangential stress is zero, but the case of a solid 
boundary is by for the most common and it alone need is 
examined in detail here. An irrotational motion of the 
fluid is determined completely by the condition of zero 
flux of mass across each element of the solid boundary, 
and the unique irrotational motion almost in-evitably has 
a non-zero tangential component of relative velocity of 
the fluid at the solid boundary. Thus, the motion that 
would be generated from rest in the absence of diffusion 
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of vorticity across the boundary of the fluid is 
accompanied by a non-zero tangential relative velocity to 
be zero at «each point of the solid boundary, however, 
small the viscosity may be, the vorticity in this flow is 
infinite at the boundary. This sheet of infinite 
vorticity at the boundary is the source from which once 
viscosity is allowed to act vorticity diffuses into the 
interior of the fluid. 

4.4. IRROTATIONAL MOTION 

If the vorticity is zero at all points in a 
region (except certain special points, called singular 
points. Where the velocity or the acceleration is 
theoretically zero or infinite) then the flow in that 
region is said to be irrotational . Flow in regions, where 
the vorticity is other than zero is said to be 

rotational. In practice there may be rotational motion in 
one part of a flow field and irrotational motion in 
another part. 

The concept of irrotational flow lies behind 

much of what follows and we may here pause to consider 

its physical interpretation irrotational flow is flow in 
which each element of the moving fluid undergoes no net 
rotation (with respect to chosen coordinate axes) from 
one instant to another. A well-known example of 

irrotational motion is that of the carriages of the big 
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(ferries) wheel in a fair-ground although each carriage 
follows a circular path as the wheel revolves/ it does 
not rotate with respect to the earth the passengers 
fortunately remain upright and continue to face in the 
same direction. 

Two examples of fluid flow are/ depicted in 
Fig. 4.2. A small element of fluid as therefore 
convenience represented by a quadrilateral with gxes AB 


and CD. 

At (a) the axis 

AB 

rotates 

clockwise 

as the 

element 

moves along/ but 

CD 

rotates 

an equal 

amount 


anticlockwise so that the net rotation is zero. Although 

undergoing distortion the element is thus not rotating; 

the flow is irrotational . At (b), however/ both axes 
rotate in the same direction/ there is thus rotation but 
little distortion. 

Rotation of a fluid particle can be caused only 
by a torgue applied by shear forces on the sides of the 
particle. Such forces/ however/ are impossible in an 
ideal. Inviscid/ fluid and so the flow of an ideal fluid 
is necessarily irrotational. 

To illustrate the difference between rotational 
and irrotational flow by a simple example. Consider some 

fluid rotating like a rigid body around on axis at right 

angles to the place of the paper with an angular velocity 
(Fig.4.3). At radius r^, the tangential velocity 
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U 2 = '^^2' radius the velocity is = wr^. The 

circulation along the boundary of the cross-batched 
region in Fig. 4.3(a) can be computed easily, since the 
two radial lines do not contribute to it. 

Therefore, 


ui r^ 0 


dr 

) 

^ = 2ui A 

In other Words, the circulation is equal to the 
product of the angular velocity and the enclosed area 
A. Obvisiously, and circulation is valid for any 
arbitrary circular sector, even an infinitesimally small 
one, since the angular velocity is the same 

everywhere. The arguments regarding the rotational vector 
show that its component must be identical with (of 
course, and are zero in this case), so that all 
particle, in the fluid have rotation and the rotational 
vector has a uniform value throughout. 


r = uj ''2 ^ 


= 0w (r2 


The area of the region is 


r ^ (rjZf) 


2 


f 2 2 

(r2 - rf 
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To contrasty consider a flow system in which the 

stream lines are also concentric circles around O but in 

which the tangential velocity is inversely proportional 

to the distance from the origin or equal to C where C is 

r 

a constant, if the circulation is computed around the 
area similar to that in Fig ,4. 3, then 

^ - “iv • 0 

i.e., the region is irrotational . It is easy to show that 
this is true for any arbitrary line abed which does 
not include the centre of rotation 0, Fig. 4. 3(b). On 
other hand; the circulation of a closed circle around O 
is 

r->=2TNr. C = 2T\C = constant 

\ 

This circulation is independent of r and gives 
physical meaning to the constant C. Since is different 
from zero, dotation must exist at point 0. Therefore, 
point O is a singular point at which u is infinitely 
large . 

Now we will consider the properties of a 
velocity field v satisfying the equations 

V.v = 0, ^x V = 0 (4.21) 

We can see that most fluids behave under a wide 
range of flow conditions, as if they were nearly 
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incompressible, and also that flow fields with the 
seemingly restrictive property of zero vorticity over 
large parts of the field area, for dynamical reasons, 
remarkably common. Study of irrotational solenoidal 
vector fields therefore has great practical claue in 
fluid mechanics. 

The implicity of the equation (4.21) has also made 
possible extensive mathematical of powerful analytical 
techniques . 

In a fluid in which the instantaneous 
distribution of velocity is v (x) material elements are 
being subjected to translation and a pure straining 
motion without change of volume and without superposed 
rotation 

Since V v is zero at all points of the fluid. 
Stoke 's theorem show that 

0 V. dx = 0 (4.22) 

for all reducible closed curves lying within the fluid, 
because it is always possible to find an open surface 
bounded by any such reducible curve and lying entirely in 
the fluid. If 0 and P are two points in a connected 
region of fluid and and are two different curves 
joining 0 and P in such a way that the two together from 
a reducible closed curve lying entirely in the fluid, 


( 97 ) 


we see from equation (4.22) that 


V 


dx 


dx 


C 


1 


C 


2 


The line integral of v over a curve joining O to P and 

lying within the fluid thus has the same value for all 

members of a set of paths of which any two make a 

reducible closed curve and depends only on the position 
vector x^ and x of O and P respectively. It is therefore 
possible to define a function 0 (x) such that 

c:^(x) = 9^(x ) + (T dx 

' ° J (4.23) 

o 


in which the integral is taken over one of the paths in 
the set mentioned. The vector gradient of (x) is found by 
varying the position of P, giving 

V (x) = V (x) (4.24) 

is termed as the velocity potential for the field v 
(although there is no question here of an interpretation 
of ^ as a potential energy function) . It is customary 
to leave the position unspecified, since the 
difference between the values of corresponding to two 
different choiceoof x^ is independent of x and so without 
effect on 

V (x) 
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The introduction of the function by means of the 
relation (4.24) ensures that the equation ^ x "v = 0 is 
satisfied identically, and the three unknown scalar 
components of v are thereby determined by a single 
unknown scalar function ^ . The first of the equation 
(4.21) then requires 





0 




(4.25) 

at all points 

of the fluid. 

This 

equation 

for 0 , 

known 

as Laplace's 

equation, appears 

in many 

branches of 

mathematical 

physics, and 

many 

general 

results 

about 

functions satisfying the 

equation are 

known 

. The 


linearity of the equation is note-worthy and accounts for 
the relative simplicity of analysis of irrotational 
solenoidal * flow, the dynamical equations governing the 
change of the velocity distribution in a fluid from one 
instant to the next are in general nonlinear but in the 
particular case of irrotational solenoidal flow the 
constraints on the velocity distribution are so strong as 
to require the spatial distribution of v to satisfy the 
simple linear equations (4.24) and (4.25) independently 
to temporal changes. 

4.5 This Section contains my two research papers. 


PAPER- I 
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ON THE VORTICITY OF MHD UNSTEADY HELE-SHAW FLOW OF 

NON-NEWTONIAN FLUID 


ABSTRACT 


In this paper an attempt has been made to study 
the vorticity of unsteady flow of non-Newtonian fluid of 
Kuvshinki's type under the influence of uniform 
transverse magnetic field, assuming the pressure gradient 
to be proportional to exp (-nt). The effect of magnetic 
field parameter (M) and relaxation time (T) on the 
vorticity of flow have been studied numerically and 
graphically . 


INTRODUCTION 

e 

The steady Hele-Shaw flows have been discussed 
by Lamp [5], Reigels [7], Thompson [10], Lee and Fung [6] 
and Buckmaster [1]. In most of these investigations, the 
pressure gradient was assumed to be consent. Swaminathan 
[9] has investigated the unsteady Hele-Shaw flow of 
viscous, incompressible fluid taking the pressure 
gradient to be a function of time. Gupta et al . [2, 3] 
studied the same problem for viscoelastic fluid. The 
purpose of this paper is to obtain an expression for 
vorticity of unsteady Hele-Shaw flow of non-Newtonian 
fluid under the influence of transverse magnetic field. 
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The effect of various parameters viz., magnetic field 
parameter (M) and relaxation time (T) on the vorticity of 
flow have been discussed numerically and graphically. 


FORMULATION AND SOLUTION OF 
THE PROBLEM 


Let us consider the flow of a non-Newtonian 

fluid confirned between two parallel plates Z = + h past 

2 2 2 

a circular cylinder x +y =a, h^Z<h under the 
influence of uniform transverse magnetic fluid. 


The non-dimensional equations governing the 
motion of non-Newtonian fluid of Kuvshinski type in the 
Hele-Shaw cell (Gupta et [2, 3] after applying the 

uniform transverse magnetic field of intensity Ho, under 
the assumption that the effect of the induced magnetic 
field and electric field produced by the motion of the 
electrically conducting fluid is negligible and no 
external electric field is applied, are given by. 






St ’St.’" 


( ^ -t-T^ 1 ^ O 


CM -W 

I M ■'^ 3 / 


^ '■A 
tD •><- 


-V 


C)T> 

2 


Where 


^ V 

I)p0o 

-7 

•7- ~ 


I iui ; 




■H 




cx 


5 Oo 

5 0 '- c. 


T 


( 


Xo"^ 


^ -4^ = X ^ 

) is the relaxation time parameter. 


the magnetic permeability. 

6 the electrical conductivity. 


The boundary conditions are, 


u=0,‘i^=0, on Z = + h 


(4.30) 


Using (4), the equation (1) and (2) yield 


( \irT ^ \l^X ) = 

■D t ^ -syi. -o ^1- 


Equation (4.28) shows that p is independent 

Therefore, p is a function of x, y and t. 

\jX~ 

^ and Vr 
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Where d* is some function of x & y. 
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of Z . 


(4.32) 


Substituting the values of u andi^from (4.32) in equation 
(4.29), we get, 


'bX.l- 


■t 




= 0 


(4.33) 


Again substituting the values of u and v in (4.26) and 
(4.27) and then integrating, we get, 


C'+"r|t^l= = 


(4.34) 
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Where g(t) is an arbitrary function of time. 

Let the pressure gradient be proportional to exp(-nt) . In 
equation (4.34) we assume that 

- Tf^^ - (H-MT)f^ = - Ae~”^ ..(4.35) 

Where N ^ I + and A is a given constant* 


Let f (t,z) = e"”^^ F (Z) (4.36) 

Solving (4.35) after using (4.36), under the corresponding 
boundary conditions, the function, f(t,z) is given as 


Where 


—n t 

f (t,z) = - ( 1 _ ^ ^ 

^1 

a2 = n (1+MT) - (M+n2^j 


(4.37) 

(4.38) 


The function 0 (x, y) can be obtained by solving (4.33) 

subject to the condition 


u. Co4 e -t- i> - O c>v\ "Vv- O, Ov- ^ ov>K=fOu . (4.39) 

c> k"” 


Where 


and 


X = r Cosd- , y = r & 




"c) W. 

Hence 0 (x, y) 


(r +^)Coi& 


(4.40) 


From equation (4.32), (4.37) and (4.40), the 

dimensionless velocity components are given by 
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The Yortlcity of flow (16) is given by 


I’C I - i-ti l 2 Lr^; r 1 


■i- ^Vn 


(4.42) 


The study the effeot of magnetic field parameter 

relaxation time (T) on the vorticity of flow, 

numerical calculations have been carried out and 

vorticity profiles have been ohi-Air>orj -p 

oeen obtained for various values 

of M and T. 


TABLE 1 

h. =1, h_=l, A = 10, M = 0, t = 0 



0.1 

0.2 

0.3 

0.4 

0.5 

1.8477 

3.6770 

5.4695 

7.2074 

8.8732 

1.5954 

3.1780 

4.7352 

6.2546 

7.7240 

1.3977 

2.7872 

4.1598 

5.5076 

6. 8254 
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TABLE 2 


ri=l , 

L 

h=l, A=10, M= 

e 

o 

II 

to 

9 

o 




T 


0.0 0.1 


jEtatm 

0.4 

0.5 

0 


0 1.3143 

2.6219 

3.9165 

5.1915 

6.4406 

0.2 


0 1.2389 

2.4730 

3.6971 

4.9065 

6.0962 

0.4 


0 1.1702 

2.3380 

3.4983 

4.6480 

5.7838 




TABLE 3 





n=l , 

j>- - 

h=l, A=10, M=0 

, t=0.2 




T 1 

Ksa 

d . 0 0.1 

6.2 

0.3 

0.4 

0.5 

0 


0 1.5126 

3.0104 

4.4781 

5.9009 

7.2648 

0.2 i 


0 1.3062 

2.6019 

3.8769 

5.12089 

6.3239 

0.4 

_j 

0 1.1444 

2.2819 

3.4061 

4.5092 

5.5856 


1 


TABLE 4 





rL=l / 

h=l, A=10, M=0 

, t=0.2 




T 


0.0 0.1 


mm 

0.4 

0.5 

0 


0.000 .1.0760 

2.1467 

3.2066 

4.2505 

5.2732 

0.2 


0.000 1.01439 

2.0247 

3.0269 

4.0171 

4.9912 

0.4 


0.000 0.9585 

1.9142 

2 .. 8641 

3.8054 

4.7354 
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DISCOSSION 

The effect of magnetic field parameter (M) and 
the relaxation time (T) on the resultant vorticity are 
shown graphically and in tabular form. From fig. 1 and 2 
and tables 1, 2, 3 and 4 it is clear that at Z=0, 

vorticity of the flow is zero. But as was shown by Gupta, 
Sisodia and Sharma (3) the velocity components at Z=0 
have maximum value. Hence we can easily conclude that a 
pure irrotation motion exists in between two planes at 
Z=0. As Z increases the vorticity also increases and at 
Z = 0.i3 it-i value becomes maximum. But as the velocity 
components become zero at Z=0 . 5 , hence we have pure 
vortex motion there. 

As relaxation time increases, the vorticity 
value is depressed and its rate of increase also comes 
down. Hence the role of relaxation time is to decrease 

the vorticity of the flow. Although the increasing 
tendency of vorticity does not change but its rate of 

increase comes down. 

ASj the magnetic field is induced, the rate of 
increase is further supressed. Although at h=0.5, we get 
pure vortex motion but amount of rotation becomes lesser 
and lesser by the induction of magnetic field and 


relaxation time. 



( 106 ) 


The increase in time also suppresses the 
vorticity . 

Thus, the induction of magnetic field, 
relaxation time and time all combindly try to decrease 
the rate of increase of vorticity and region of 
irrotationallity at the central time and the region of 
maximum vorticity along the planes can be found in the 
flow. 


•kifk'k'k'k'k'k* 
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PAPER - II 

4.6 ON THE VORTICITY OF FLOW OF AN ELASTICO-VISCOOS 
FLOID NEAR AN ACCELERATED PLATE 

ABSTRACT 

In this paper an attempt has been made to study 
the vorticity of flow of an elastico-viscous fluid in 
case of Walters liquid ^ past a uniformly accelerated 
plate. The effect of elasticity parameter K on the 
vorticity of flow have been studied numerically and 
graphically. 


INTRODUCTION 

The study of the flow behaviour of Non-Newtonian 
fluid is receiving the attention of research workers in 
technological field. It has manu practical applications. 
The flow of an elastico-viscous fluid with short memory, 
past and impulsively started plate was studied by 
soundalgekar [8], Haldavnekar [4] has studied the flow of 
an electro-viscous fluid (waiters liquid fi) near an 
accelerated plate. In this paper we have obtained an 
expression for the vorticity of flow of an 
elastico-viscous fluid near an accelerated plate. Effect 
of elasticity parameter on vorticity of flow have been 
studied numerically and graphically. 


FORMDIATION AND SOLUTION OF THE PROBLEM 


Let us suppose that an infinite plate is 
uniformaly accelerated in the X-axis which is taKen along 
the plate. The Y-axis is taken normal to the plate. The 

fluid and the plate are assumed to be stationary at t^ <0, 

. . / 
where t' is the time. The plate is accelerated at t=0. 

The equation governing the motion is 
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The 

boundary conditions are 


\ 

u ■ 

= 0 

at t' 0 

1 

1 


u = 


at y' = 0 

{ 

1 


! 

U = 

0 

at y ' 




For 

such fluids Kjj 

<<1 

t 



To solve the equation (i), let us assume that 
u = u^ + K ' u, . . . . 


(4.43) 


(4.44) 


(4.45) 


putting the value of u' from (4.45) in equation (4.43) & 

(4.44) and equating the coefficients of different powers 
of K' neglecting those of , we get 
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u„ = 

Uj^=0 at y'=0 1 

Uo = 0 

u^=0 at y'-^o:? J 


(4.47) 


(4.48) 


Applying Laplace transform technique, the solution of 
equation (4.46) & (4.47) was given by [4] as, ^ 
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Introducting the non-dimensional quantities. 


-rre' 


A^' 


we have from equation (4.45) adn (4.49), 




VL ^ ' + '^) -t 


(4.50) 


The vorticity of flow (4.50) is given by 

-v- S r.tl ....( 4 . 51 ) 
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The study the effect of elasticity parameter K 
on the vorticity of flow (4.51) numerical calculations 
have been carried out and vorticity profiles have been 
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obtained for various values of K 


K 

ESI 

o 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 

0 

i 


'-0.5642 

0.2834 

1.0774 

1.6157 

1.7703 

1.5685 

1.1697 

0.5 1 

1 


-0.7052 

+0.1675 

1.02251 

1.6257 

1.8222 

1.6313 

1.2163 

1.0 


-0.8463 

1 

0.0449 

1.0774 

1.6358 

1.8741 

1.6942 

1.2683 


CONCLUSIONS : 

The variation of vorticity, and the influence of 
elastic parameter K on it has been depicted in the 
firgure 3. From the vorticity profile and the numerical 
calculations, we observe that 

i 

(i) The vorticity increase with ■v\_and somewhere near 
the plate (0.5 >'^> 0) region of zero vorticity 
is observed. 

(ii) The vorticity increase upto "'^=2 and then it 
begins to decrease. 

(iii) As the value of K increases, the value of 
vorticity also increases. 


★ ★ Tfc- * ★ ★ ★ 
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CHAPTER V 


CHAPTER - 5 


VORICITY IN POROUS MEDIUM 

5 • In this chapter fluid motion in porous medium 
has been considered. The effect of permeability of 
porous medium, time and electrical conductivity has 
been studied. This chapter contain two research papers t 
I* On the Vorticity of flow of a viscous fluid in 

a porous medium due to an oscillatiny flat 
plate. 

2. On the Vorticity of Hydromegnet i c Visco-elastic 

fluid near oscillating porous flat plate. 

5.2 MHD FLOW IN POROUS MEDIA 

The study of physics of flow through porous 
media has become the basis for many scientific and 
engineering applications. The study of MHD flows through 
porous media is receiving considerable attention of 
researchers owing to its importance in many scientific 
and engineering fields, particularly in power 
generation. MHD pumps are already in use in chemical 
engineering technology for pumping electrically 
conducting fluids in the Atomic Energy centres. Such 
type of flows veimportant in petroleum technology to 
study the movement of gas, oil and water through the 
oil reservoirs and in chemical engineering for 
filteration and purification processes. Beyond this, the 
study is widely applicable in soil mechanics, water 
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purification, ceramic engineering and powder metallurgy. 

The study of MHD flow through porous media is 
also of great interest in geophysical to study the 
interaction of the geomagnetic field with the fluid in 
the geothermal region. Water in the geothermal region is 
an electrically conducting liquid because of high 
temperature. With the fuel crisis deepening all over the 
developed world, attent ion is turning to the 
utilization of the enormous power beneath the earth's 
crust in the geothermal region. 

Another potential geophysical application may 
be in the exploration of geopressured reservoirs. In 
these reservoirs, water at elevated temperature exists 
at enormously high pressure because of the weight of the 
overlying rock and the geomagnetic field. The upflowing 
water from geopressured wells can run hydraulic turbines 
to produce electricity, while the heat in the water can 
simultaneously be extracted to run steam turbines, again 
producing electricity. 

Ahmadi and Manvi [1] derived a general equation 
of motion of fluid flow through porous media and 
applied to some basic flow problems. Gulab Ram and 
Mishra [5] studied MHD flow of conducting fluid through 
porous media. Varshney [30] studied unsteady MHD flow 
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through a porous medium in a circular pipe. Gupta [6] 
studied unsteady MHD flow through porous media in a 
channel of cross-section circular sector while Kumar, 
Gupta and Varshney [7] studied MHD channel flow through 
a porous medium. 

In the present thesis I have also tried to 
study two MHD flows through porous medium and this work 
has been included as two research papers. 

5,3 VORICITY IN' MAGNETOHYDRODYNAMIC FLOW 

The workers dealing with MHD have paid less 
attention on the property of vorticity of the flows. 
Recently, Mittal [11,12,13,14,15] and Mittal et al 
[16,17, 19,20, 21,22,23] have studied the vorticity 
of some MHD flows. Here, in this section, the voriticity 
of MHD flow on or around a vertical porous plate under 
different magnetic field has been studied. The effect of 
the applied magnetic field on the vorticity of the flow 
has been shown by 'vorticity profiles' under different 
contions. Attempts have been made to find the conditions 
of irrotationality of the motion which makes the study 
of flow simple. Also, in the region of irrotationality 
the energy losses minimize. 

Phenomenon of sueprposability and self 
superposabil ity having been proved to be the most 
applicable techniques in solving the hydrodynamic and 
magnetohydrodynamic problems have been frequently used 
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here in bringing about some very important conclusions 
about the MHD duct flow. 


In the following pages/ I have given two 
research papers/ which are my own contributions and deal 
with the study of vorticity of different MHD flows under 
different conditions. 


5.4 ON THE 

VORTICITY 

OF 

FLOW OF A VISCOUS 

FLUID IN A 

POROUS 

MEDIUM DUE 

TO 

AN OSCILLATING FLAT 

PLATE 


ABSTRACT 

In this paper an attempt has been made to study 
the vorticity of unsteady flow of a viscous fluid in a 
porous medium due to an oscillating flat plate. The 
effect of permeaility of the porous medium and time on 
the vorticity of flow have been studied numerically and 
graphically by plotting the vorticity profiles and 
results have been discussed critically. 

INTRODUCTION 


Fluid 

importance 
technology , 
hydrology/ 
industrial 


flow through porous media is 
in a wide range of disciplines 
6.g./ in soil mechanics/ 
petroleum engineering/ water 
filtration/ ceramic engineering 


of considerable 
in science and 
ground water 
puri f icat ion , 
etc. 


media 
Da rcy 


The 

theory 

of 

laminar fl 

ow through/ 

homogeneous 

i s 

based on 

an 

experiment 

originally 

conducted by 

[3] 

/ Muskat 

[5] 

discussed 

the flow til 

irough porous 
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media in connection with filtration. Ahmadi & Manvi [1] 
gave a general equation governing the motion of a viscous 
fluid through porous media. The flow of a viscous fluid 
about an infinite flat wall executing harmonic 
oscillations parallel to itself was first studied by 
stokes [11] • Rukmangadachari [10] has obtained an 
expression for velocity of flow of a viscous fluid in a 
porous medium due to an oscillating flat plate. 

In this paper we proposed to study the vorticity 
of the flow of an incompressible viscous fluid through a 
porous medium due to the oscillations of a flat plate in 
its own plane. Numberical calculations have been carried 
out and vorticity profiles have been obtained for values 
of k, the permeability of the porous medium and time. 

FORMULATION AND SOLUTION OF THE PROBLEM 

The equation governing the motion of an 
incompressible viscous fluid permeating a porous medium 
of premeability k by [1] are 

V-l? - 0 (5.1) 

= - 'Vp-^ + ^ it (5.2) 

where 1? is the velocity vector/ p is the density / is 
the coefficient of viscosity and p is the pressure. 

Let us consider an infinite flat osci 1 lat ing . in 

i , ^TT t . ... 

its own plane with velocity u^e * in an infinite 
expanse expanse of an incompressible viscous fluid 
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permeating a porous medium of permeability k. Take the 
X-axis along the plate and the y-axis normal to the plate 
so that the plate is defined by y=0. For rectilinear 
flow, the non-zero component of velocity, u is a function 
of y and t only. The equation (5.1) is identically 


satisfied and equation (5.2) becomes 


P __ -V 


C T.) 


Xvx e-ooun- Uv - ‘a ^ ) e 
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In real terms the velocity is given by 

^ «r^ Cc^> C<rt 5 ^ ^ 

\, 
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Where 


a 



(5.9) 


i ^ I 


The vorticity of flow (5.8) is given by 

“5 - ^ __ 


Ua- _ 
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UoTfi, ' ^ ^ 

The equation (5.10) reduces to 


^ \ -V C \-V ^ ^ ^ * 


We have 


(5.10) 


: . ( 5 . 1| 
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TABLE 1 

t = 0 . 4 ^ 


m 


0.0 

o '. 2 

6.4 

0.6 

0.8 

1.0 

1.2 

0.2 

< 

0.9944 

0.8134 

0.6653 

0.5442 

0.4451 

0.3641 

0.2978 


> 








0.6 


0.9618 

0.7813 

0.6346 

0.5156 

0.4188 

0.3402 

0.2758 

1.0 


0.9211 

0.7399 

0.5943 

0.4863 

0.3835 

0.3081 

0 . 2-474 


TABLE 2 


t = 0.8 



IPj 

0.0 

0.2 

0.4 

0.6 

00 

• 

o 

o 

« 

1 — 1 

1.2 

0.2 


0.9938 

0.8127 

0.6648 

0.5437 

0.4447 

0.3638 

0.2975 

0.6 


0.9599 

0.7797 

0.6334 

0.5145 

0.4179 

0.3395 

0.2758 

1.0 


0.9184 

0.7377 

j 

0.5926 

0.4849 

0.3824 

0.3071 

0.2467 


TABLE 3 


t = 1.2 


k 

3 _ 

o 

o 

0.2 

0.4 

0.6 

00 

o 

1.0 

1.2 

0.2 


0.9928 

0.8121 

0.6643 

0.5433 

0.4444 

0.3635 

0.2973 

0.6 


0.9579 

0.7781 , 

0.6321 

0.5135 

0.4171 

0.3388 

0.2752 

1.0 

-1 

0.9156 

0.7355 

0.5908 

0.4834 

0.3812 

0.3062 

0.2460 
































( 125 ) 


CONCLUSION 


The effect of permeability on the vorticity of 
flow of a viscous fluid in a porous medium due to an 
oscillating flat plate has been studied through the 

calculations given in tables 1-3 and figure 1-3. 

It is clear from above observations that the 
vorticity is maximum at the plate which is 

self-explanatory because of oscillations in it. As we 
move away the plate the value of vorticity lessens but we 
may practically not expect any region of irrotationality 
in its vicinity. With increase in time although the 

vorticity lessens but very slightly. Also the increase in 
permeability of the porous medium does not have any 
marked effect on it. Although, with the increase in 

permeability the net value of vorticity decreases 

slightly, but its mode of decrease with y* remains almost 
same. So why, the vorticity profiles are almost parallel. 

5.5 PAPER - II 

ON THE VORTICITY OF HYDROMAGNETIC 
VISCO-ELASTIC FLUID NEAR OSCILLATING POROUS FLAT PLATE 

i 

ABSTRACT 

In this paper an attempt has been made to find 
the vorticity of Hydromagnetic visco-elastic fluid 
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( Kuvshinki ' s fluid) with small electrical conductivity 
near an infinite insulated porous flat plate oscillating 
hormonically in its own plane in the presence of a 
transverse magnetic field of uniform strength fixed 
relative to the fluid. Laplace transform technique has 
been used for the solution of the problem. Small uniform 
suction has been imposed along the surface of the plate. 


Numerical 

calculations 

has been 

carried 

out 

for the 

vorticity 

field for 

different 

values 

of 

suction 


parameter, magnetic parameter and elastic parameter at 
T=0 and T= ('/2 and vorticity profiles obtained and result 
discussed critically. 

INTRODUCTION 

Stokes [12] and Rayleigh [8] investigated the 
flow about an infinite flat wall executing linear 
harmonic oscillations parallel to itself. Rossow [9] 
discussed the impulsive motion of the infinite plate in a 
viscous incompressible conducting fluid in the presence 
of an external magnetic field. Ong and Nicholls [6] 
extended the method to obtain the flow of an electrically 
conducting fluid near an infinite flat wall which 
oscillates in its own plane. Chaudhary [2] has given a 
detailed study of the motion of a viscous incompressible 
and electrically conducting fluid near an infinite 
oscillating flat plate in the presence of an uniform 
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magnetic field relative to the fluid. 

In this paper we have made an effort to find the 
vorticity of the hydromagnetic flow of a visco-elastic 
fluid near an oscillating porous flat plate. The effect 
of various parameters on the vorticity of hydromagnetic 
flow has been discussed numerically and graphically. 

FORMULATION AND SOLUTION OF THE PROBLEM 

Let us consider an infinite insulated porous 
flat plate executing linear harmonic oscillations in its 
own plane in presence of an externally applied transverse 
magnetic field of uniform strength fixed relative to 

the fluid. As the plate is infinite in length and uniform 
suction is imposed over it, the physical variables depend 
only on y, the coordinate perpendicular to the wall and 
t, the time. The pressure p in the fluid assumed 
constant . 

(5.12) 

(5.13) 




The equations governing the flow are 








X , X. 

e U. 




^ o (Continuity equation] 


From (5.12) and (5.13) we have 


(5.14) 
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Her 0 is the relaxation time, signifies that the shear 
stress decays are exp (“t/*^ ) when the motion stops 
suddenly . 


If V„ represents the velocity of suction or injection at 
the plate, the equation (5.13) with the condition y=0, 
V=Vo leads to the result V=Vo every where. 


Now eleminating , between 


(5.12) and (5.13), We get 


[l 4-> < ^ • 15 ) 

Where u is the velocity parallel to the plate, p the 
density, ^ the kinematic viscosity, o- the 

electrical conductivity and hi the magnetic permeability. 

The term l.'VvWA. •= ) represents the Ponder 

motive force parallel to the plate. Here the magnetic 
Reynolds number R s' is considered small so that the 
induced magnetic field due to the flow may be neglected 
with respect to the applied magnetic field. The boundary 
conditions of the problem are 


n- 


<j tk — ^ Co^ 'vv^ 


Ciy 


(5.16) 


Let us assume that 


M 








'si -i:! ^r- . ‘\J lo 


(5.17) 
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Accordingly the equation (5.15) reduces to 




•ST 




^•vv 


'X 


C)T 


-vV. 


■a-vy 




Sv- 


(5.18) 


The transformed boundary conditions are : 


'VV^ 0 


VA — > 6 


(5.19) 


Initially, it is assumed that, the fluid is at rest, 
i .e. , 

•V\.-o C) tvT’ Tr-O 

^ -ST 

Applying Laplace Transform on each term of equation 
(5.18) we get 

- cXG: lV + M+V"-Y\'rMV2S == o 

Aj> -- - ; . 


(5.20) 


Where 


Va 


°^u, C-Vt) AT 


The solution of equation (5.20) is given by 


Jb} -S-Y^L" ‘^■^r^CV+j?>)( WV^"] (5.21) 


Y Vo^ 41 , 7 1 


j"; 


(5.22) 
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Using Laplace inversion, theorem, we get 

lx. ,T ) = ^ ( __L ^ ^ A , 

T-ttX \+v^ 

(5.23) 


Using the theory of residues, the velocity distribution 
\A(^X,T) is given by [4] as. 


-'^o'^^V\rVV^ S"-VC^ Cj6 I'b -V <1. P StMft))..( 5.24) 


Where V^. 

f Chp<J> 6 - 

0 ^ ^ 

The vorticity of flow (5.23) is given by 

^ =• ^ -S--^V \ ^ ^ S' -V cv f 

^(> S -V G. f Lc^ ®)Ca5 (5^25) 


CASE-I When vy=3.e> ^ ^Tv^ri-D 

We have fi ;> ^ 

The velocity field (5.24) reduces to 


..(5.26) 
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and the vorticity of flow (5,26) in this case becomes 


0>t 








( 5 - 27 ) 


CASE-II When 


V^C) ^ V - O ^ n 0 

ijh ^ CO ^ Y 

and the velocity field in this case is 

(5.28) 


y- 

- 


The vorticity of flow in this case is given by 

® V ^ \ 


ijSow.'^) V M ) ^ \t- 'iy ttW) (5. 


29 ) 


CASE-III When 5 V, ^ 'f-o.o ^ 

p. -> M 


we have 


and 


X 
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and the velocity field in this case is given by 

^ lx. - u 


-X. ( T' \ 




:{5.30) 


The vorticity of flow in this case is given by 
< c r; rr: 

'^o 2. 

o 


X 

-V4H 




<Lto/s t,T '-IVo 




X (t->= 


Vsi-' 

— t~. , r^L 


<° + C J 


- {.Mo j 


(5.31) 


't = 0, M = 0, T = 0 


TABLE 1 


"V^g^nO O TTo TTs JTo ITS TTo" 



0 

2 


1 


-1.0000 -0.9018 -0.6952 -0.4712 -0.2431 -0.2008 -0.1648 
-0.1408 -0.1720 -0.1517 -0.1161 -0.0814 -0.0535 -0.0333 
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TABLE 2 

^'-0, M-0, T- 7r/2 



@9 

0.0 


1.0 

1 .5 

2.0 

2.5 

3.0 

0 

B9 

0.000 


0.3203 

0.3021 

0.2401 

0.1674 

o . io : 

2 

E9 

! 0.645 

0.2445 

0.1911 

0.0923 

0.0387 

0.0117 

- 0 . 00 ( 


f-0, M-2, T-0 Table 3 




'EBIH 

mm 

1.0 

1.5 

2.0 

2.5 

3.0 

0 

WBM 

- 2.0582 

- 1.0196 

- 0.4903 

- 0.2288 

- 0.1034 

- 0.0450 

- 0.0187 


B1 

DB 

- 0.7639 

- 0.5348 

- 0.3668 

- 0.2465 

- 0.1709 

- 0.1043 

1 


TABLE 4 


f»0, M-2, T- Tr /2 




0.0 

0.5 

1 .0 

1.5 

» 2.0 

2.5 

3.0 

B 

m 


0.0613 

- 0.0550 

- 0.0726 

- 0.0506 

- 0.0270 

- 0.0192 


m 

0.4028 

0.1694 

0.0406 

- 0.0245 

- 0.0523 

- 0.0595 

- 0.0561 


TABLE 5 


F - 2 , M-0 , T-C 
1 


E 


0.0 

0.5 

1 .0 

1.5 

2.0 

2.5 

3.0 


m 


2.0265 

1.7087 

- 0.2418 

- 1.4964 

- 1.2095 

- 0.012 9 


m 

- 0.5068 

- 2.3169 

- 2.6871 

- 1.7011 

- 0.0541 

1.1629 

1.5 508 


TABLE 6 




0.0 

* U 5 i ® ' V 

0.5 

1.0 

1.5 

2.0 

2.5 

3.0 



- 4.1163 

IHH^I 

0.9067 

1.6709 

0.5708 

- 0.7830 

- 0.8301 


mssM 

- 3.0911 

- 1.7869 

0.2066 

\ 

1.7488 

1.8751 

1.2565 

0.2 951 
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TABLE 7 


V -2, M^-2, T-0 




0.0 

0.5 

1 .0 


2.0 

2.5 

3.0 

0 

WSi 

-4.4721 

0.6751 

0.6711 

0.1857 

0.0617 

-0.0235 

-0.0027 

2 



0.7777 

0.0900 

0.0173 


-0.0540 

-0.0315 


TABLE 8 


'f-2, M^ = 2, T= -Tr/2 


% 


1 

0.0 

0.5 

1 .0 

1.5 

2.0 

2.5 

3.0 

0 


-4.4721 

-0.2478 

0.5779 

0.0779 

0.0217 

0.0195 

-0.0005 

2 

n 

-3.0673 

-1 •8464 

-0.2129 

-0.0214 

0.0141 

-0.0139 

-0.0265 


The effect of suction, magnetic field parameter and 
elastic parameter on the vorticity of the flow, at T«0 and 
T“ '^12 has been shown in figures. It is clear that 
variation of vorticity is not uniform in the flow and it 
markebly depend on ^ M and T. In almost all the cases 
the vorticity tends to become zero. In some cases the region 
of irrotationality lies near the plate (the distance 
depends on the parameters) and in others it lessens as we 
move away from the plate. Obviously the elastic parameter 
and the suction or injection velocity govern the vorticity 
of flow while the introduction of magnetic field adds more 
complications to it. 
















































Vo =o. M = q,T=^? 


Lj 










^Vc. = o. y = 2, M=0, r^'^Ai 

.Vo = 2,f-=5,M-o. T-^Va 


$ ( Toi4 S’- 




e 
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